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Abstract 


A two-dimensional to three-dimensional global/local finite element approach was developed, veri- 
fied, and applied to a laminated composite plate of finite width and length containing a central 
circular hole. The resulting stress fields for axial compression loads were examined in detail for 
several symmetric stacking sequences and hole sizes. Verification was based on comparison of the 
displacements and the stress fields with those accepted trends from previous free edge investigations 
and a complete three-dimensional finite clement solution of the plate. Hole diameters of one, three, 
and six inches in plates 18 inches long, 12 inches wide, and 0.1 inches thick were considered. The 
laminates in the compression study included symmetric cross-plv, anglc-plv and quasi-isotropic 
stacking sequences. The entire plate was selected as the global model and analyzed with two- 
dimensional finite elements. Displacements along a region identified as the global/local interface 
were applied in a kinematically consistent fashion to independent three-dimensional local models. 
Local areas of interest in the plate included a portion of the straight free edge near the hole, and 
the immediate area around the hole. 

It was found that the global/local interface should not be placed inside or through any region where 
the stress field exhibits three-dimensional effects. Interlaminar stress results obtained from the 
global/local analyses compared well with previously reported trends, and some new conclusions 
about interlaminar stress fields in plates with different laminate orientations and hole sizes are pre- 
sented for compressive loading. ILie effectiveness of the global/local procedure in reducing the 
computational effort required to solve these problems is clearly demonstrated through examination 
of the computer time required to formulate and solve the linear, static system of equations which 
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result for the global and local analyses to those required for a complete three-dimensional formu- 
lation for a cross-ply laminate. 

The Testbed, which is under continuing development by the Computational Structural Mechanics 
Group, now the Computational Mechanics Branch, was used throughout this investigation. 
Specific processors used during the analyses arc described in general terms herein The application 
of this global/local technique is not limited to this software system, and was developed and de- 
scribed in as general a manner as possible. The methodology developed is thus applicable to other 
large-scale structural analysis systems. 
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1.0 Introduction 


Since their initial use by the aerospace industry, modem composite materials have been utilized in 
a wide range of applications. Extensive testing was done initially to establish the reliability and 
service limits of specific parts. As understanding and confidence in composite materials grew, it 
was desired to eliminate some of the time lag and expense of fabrication and testing by using ana- 
lytical techniques. The inherent complexity of composite materials has meant that computer- 
implemented analysis techniques are required to predict the structural response of even simple 
components. Composite material constitutive relations are often highly coupled between extension, 
shearing, and bending, and not necessarily constant throughout a structure. Geometric disconti- 
nuities, large stress gradients, large displacements and rotations, and material nonlinearities may also 
be present. Different means of analysis exist for response prediction with varying levels of accuracy, 
capability, and efficiency. Available approaches include analytical, variational, finite difference 
methods, finite element methods, and boundary element methods. Each approach has advantages 
and disadvantages. Some of the more computationally expensive techniques have very broad 
modeling capabilities, while other less complex techniques work well only for specific problems 
with simple geometries. The analyst is required to select, or possibly develop, an analysis strategy 
which best suits the problem at hand. 
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To analyze general problems, many have turned to the finite element method. The initial approach 
adopted was to use the computationally least demanding means possible: two-dimensional ele- 
ments. For composite laminates of finite size subjected to external loads, application of a two- 
dimensional finite method can yield accurate results at locations away from traction-free edges and 
discontinuities where the stress states are more complex. 

Stress fields in composite laminates, in general, consist of six components which can be divided into 
two groups: intralaminar stresses (<r x , <x y , T*y) and interlaminar stresses (<Ti, t«, t^). Intralaminar 
stresses are the in-plane components, while interlaminar stresses arise during the transition from the 
internal in-plane stress state to that near a traction-free edge or discontinuity. The character of the 
stress distribution is dependent not only on the type of loading and the geometry involved, but also 
on the Poisson's ratio mismatch or shear coupling that may exist between adjacent material layers. 
All six components of stress can have large gradients during the transition from the internal stress 
state to the traction-free edge condition (Renieri and Herakovich, 1976), (Wang and Crossman, 
1977). Two-dimensional elements have been applied extensively during the study of a straight free 
edge of a laminated plate subjected to uniaxial tension or compression (Salamon, 1980). However, 
near a circular hole the transition from the internal stress state to the traction-free edge of the hole 
is considerably more complicated, and a two-dimensional finite element solution cannot necessarily 
capture this effect (Dana, 1974), (Salamon, 1980). As limitations of two-dimensional analyses be- 
came known and more powerful computers became available, the use of three-dimensional finite 
elements for detailed stress analysis of composite structures increased (Rybicki, 1971), (Dana, 1974), 
(Dana and Barker, 1974), (Barker et al, 1974), (Renieri and Herakovich, 1976), (Stanton et al, 
1977), (Griffin, 1983), (Griffin, 1989), (Thompson and Griffin, 1990), (Thompson et al, 1990). For 
many applications, however, a full three-dimensional analysis may be inefficient or even intractable. 
It is logical to streamline the analysis so that adequate response predictions are obtained with 
minimum computational and analyst effort. 

A simplified analysis on the entire body, or global region, is done here with two-dimensional finite 
elements. A more detailed and independent analysis is then performed on a local region of interest 
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with three-dimensional finite elements. In particular, this two-dimensional to three-dimensional 
global/local finite element analysis technique is applied to a finite laminated composite plate with 
a central circular hole, with the local regions of interest being the area around the hole and a portion 
of the straight free edge near the hole. The interlaminar stress state in these regions is examined in 
detail for several stacking sequences and three different hole sizes under compressive loading. 
Cross-ply, angle-ply and quasi-isotropic laminates in compression are considered herein. All plates 
had dimensions of 18 inches in length, 12 inches in width, a thickness of 0.1 inch, with one-, three-, 
and six-inch diameter holes centrally located. The results obtained are compared to known infor- 
mation about interlaminar stress states and to complete three-dimensional finite element solutions, 
whenever available. The development and verification of global/local finite element analysis was 
considered to be the primary focus of this particular work, and the interlaminar stress states were 
examined during the verification of this technique. In addition, the effectiveness of this computa- 
tionally more efficient method of analysis was demonstrated. 
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2.0 Literature Review - Global/Local Analysis 


Although the concept of global/local analysis is not new, little work has been published which ex- 
plains the development and verification of specific techniques, especially as applied to stress analysis 
of composite structures. Recently, more analysts have been addressing the issue of reducing the 
computational effort associated with nonlinear and transient structural analysis. It is believed that 
a considerable amount of global/local analysis has been utilized in the aircraft industry. The ma- 
jority of this information, however, is proprietary and does not appear in the open literature. This 
fact complicates the transfer of technical information about global/local analysis techniques and 
hinders careful development and verification of the methods. 

The definition of global/local analysis varies from one analyst to another, and a wide range of such 
analyses can be performed in the analysis and design of various structural components of an air- 
plane. For example, the fuselage, skin, stiffeners, and fasteners would require different levels of 
modeling detail and analysis. Concurrent work on several areas of a complex problem could occur 
more readily with global/local analyses. However, before any global/local analysis technique is 
used, development and verification of that method must be complete. Reddy (1985) gives several 
candidate problem areas that lend themselves to global/local analysis and may be used in the de- 
velopment of a particular global/local analysis technique. These include free-edge stress and contact 
stress problems, impact, fracture mechanics, and unbounded domains. All of these problems share 


Literature Review - Global/Local Analysis 


4 



some of the following features: stress concentrations which require detailed modeling, often three- 
dimensional in character; large rotations/strains which require a nonlinear analysis; local disconti- 
nuities, either geometric or material in nature; and possible material nonlinearities. Some of these 
may affect a small section of the structure in question, and attempting to incorporate all of the de- 
tails in a single model of the entire domain can lead to large, prohibitively expensive models. It is 
proposed to perform an accurate global model which determines the local effect on a gross scale, 
and then subsequently perform as many local analyses as required. These local analyses would 
account for the local effects either by further refinement of the global theory or by applying a new 
theory in the local region. Reddy (1985) and Noor (1986) identify several problem areas that need 
to be addressed for the development and use of global/local analyses. These include the selection 
of the regions and the method to be used in those regions, interfacing between the regions and the 
possibly different methods while maintaining displacement compatibility and traction equilibrium, 
and how to best compensate for local regions that are transient and may change in geometry or the 
required analytical technique during the course of the analysis. 

Various combinations of methods have been proposed for global/local techniques and for both the 
global and local areas, including analytical, variational, finite difference, finite element, and bound- 
ary elements (Reddy, 1985). The combination of different methods has been termed hybrid 
global/local analysis methods. Substructuring is the reduction of a complex structure to smaller 
parts which will allow for more detailed analysis of the localized effect. This substructuring can 
involve the use of different methods between the regions or merely a refinement of the smaller re- 
gions. Zooming is a type of substructuring when a more refined local model is generated from a 
coarse global model. As can readily be seen, global/local analysis techniques can be created which 
fall into more than one category. A comprehensive literature review of global/local analyses can 
be found in Ransom (1989). Discussion of some of the techniques applicable to structural analysis 
of laminated composites with holes will be presented herein. 

Bradford et al (1976) and Dong (1983,1985) have developed a hybrid technique, which they call 
Global Local Finite Element Method (GLFEM). This particular method focuses on sets of 
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problem-dependent analytical functions to model the global response. The resulting displacement 
and rotation field is applied to a local two-dimensional finite element model. Several problems have 
been solved using GLFEM with success: the axisymmetric scattering of elastic waves, determining 
the frequency of simply-supported, rectangular, isotropic plates with and without holes, and crack 
tip analysis in fracture mechanics. However, this method is highly dependent on the global func- 
tions and requires that the selected functions model the global behavior accurately. In addition, the 
solutions of these global and local models are coupled, requiring a reanalysis for each problem pa- 
rameter change. This method lends itself well to problems of regular geometry, but not to problems 
whose geometry or response cannot be described readily in a functional manner, such as composite 
laminates with discontinuities. 

The linear constraint method is attributed to Schwartz (1981) and was used by Kelley (1982) to 
impose the displacement field from a coarse structural model on a refined local model of a discon- 
tinuity. To apply displacements at nodes which are not in the global model, interpolation of the 
global displacements from nearby nodes was done with the global element shape functions, which 
were linear. Depending on the refinement level, the linear constraints on the additional local nodes 
could over-constrain the model and introduce errors into the displacement field. These small errors 
in the displacements could result in more substantial errors in the stress field; particularly, if stress 
gradients were present. Another problem with this method is that re-analyses are required for 
changes in the local area of interest, and that foreknowledge of the extent of the localized area is 
required prior to the global solution. 

Hirai and colleagues (1984) have developed 'an efficient zooming technique.' A coarse finite ele- 
ment mesh was analyzed, and zones of interest were identified and subdivided further. The dis- 
placements from the coarse mesh are applied as input to the finer mesh, but the effect of the refined 
region is ignored in the overall response of the coarse model. Static condensation is then performed, 
and exact structural reanalysis methods are used to analyze the new system. Successive application 
of this technique requires the information from all of the previous levels of zooming in order to 
generate the solution, and thus this method may become rather unwieldy. Further developments 
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led to the 'exact zooming method' (Hirai, 1985). The improvements eliminated some of the infor- 
mation required for each zoom level since successive zooms involved the previous level. The ex- 
panded stiffness matrix contains information from all of the previous levels in the stiffness of the 
boundary of the local area, and is used to solve only for the information at the desired location 
contained within the final refinement. Only an isotropic plate with a hole was considered. 

Global/local modeling was done by Knight and colleagues (1989) to calculate the detailed stresses 
around discontinuities in a buckled composite component with analytical failure prediction tech- 
niques. Two-dimensional finite elements were used for a coarse global model, and the displace- 
ments and rotations at the global/local interface were applied to a refined two-dimensional model 
created with a zooming technique. It was noted that a three-dimensional modeling capability would 
be required in order to determine through-the-thickness state of stress for this particular case. 
Four-node quadrilateral elements were used throughout the modeling process. 

Ransom et al (1989a, 1989b) used a two-dimensional global model to model the overall behavior 
of an isotropic panel with finite-width effects, and a composite blade-stiffened panel with a dis- 
continuous stiffener. The local models are two-dimensional, but are more refined to determine the 
local detailed stress states for specific regions. Application of the displacements and rotations to the 
local model was done with interpolation splines. The coefficients of the splines were determined 
from the global-displacement field, and so this technique does not require a priori knowledge of the 
location of the local region. In this approach, those nodes in the local model which are not in the 
global model will have the displacements determined from the appropriate interpolated displace- 
ment field. The out-of-plane strain and stress components were not presented. 

The specified boundary stiffness/force method (Jara- Almonte and Knight, 1988) uses the stiffnesses 
and forces from the global model at the nodes of the local model. This approach is like substruc- 
turing, but discards the information from the old region. Stiffnesses and forces represent the effect 
of the rest of the structure on the subregion through the incorporation of the global stiffness terms 
in the local stiffness matrix of the subregion model, and the application of the forces as external 
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loads on the local area. Global/local analysis was performed, and it was found that for the problem 
considered, their method was more accurate than specifying displacements, and was an improve- 
ment over the linear constraint and the zooming techniques developed by others. 

The development of a two-dimensional to three-dimensional global/local finite element analysis 
was reported by Vidussoni (1988). A global two-dimensional finite element analysis using 4-node 
quadrilateral elements was followed by a three-dimensional finite element analysis using 20-node 
displacement-based brick elements. Interfacing between the models was done by applying in-plane 
displacement from the global model to the independent local models. The effect of the presence 
of a circular hole on the straight free-edge stresses, and of the stress states around holes was con- 
sidered. A limitation of this work was the version of the software system and the available com- 
putational power. Part of the current work repeats some of the analyses performed by Vidussoni, 
using quadratic elements and extends the technique to other laminates. Studies done by Vidussoni 
for acceptable element aspect ratio limits, the number of elements per material layer for a given el- 
ement integration scheme, and the straight edge nodal distribution were used in the present work. 

A natural extension of two-dimensional to three-dimensional analysis is to use transition elements 
between the two models. Surana (1980) presented information about a family of transition finite 
elements between two-dimensional shell elements and three-dimensional isoparametric solid ele- 
ments. Linear elastic isotropic applications were initially considered, with further work being 
published by Surana (1982) which included a geometrically-nonlinear formulation for the family 
of elements. Liao et al (1988) investigated both linear and geometrically-nonlinear static and tran- 
sient analysis of laminated anisotropic shell-type structures. Comparison was done between models 
built of degenerated shell elements, three-dimensional continuum elements, and the transition ele- 
ments. Gong (1988) also derived a family of transition elements of variable number of nodes ena- 
bling three-dimensional finite elements to be joined to two-dimensional finite elements. 
Applications examined included an aircraft hatch door with a polymethyl-methacrylate shell. Be- 
cause of the formulation of the transition elements, their use is appropriate only in regions where 
out-of-plane stresses, (a r , r r G ) are negligible, and gradients are relatively low. 
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Several global/local techniques have been developed specifically to reduce the time and effort asso- 
ciated with solving nonlinear problems. Clough and Wilson (1979) reported work in applying a 
global/local substructuring approach to a structure with a materially nonlinear response region. 
Most of the structure was assumed to be linear with only small areas of nonlinearity. Subdivision 
of the structure into independent linear and nonlinear areas was based on the assumption that each 
nonlinear section was also independent of other nonlinear sections. During each step of the analysis 
the nonlinear subsystem stiffnesses were updated, while the evaluation of the linear stiffness oc- 
curred at the first load step, and then these were assembled to determine the total system stiffness. 
Specific examples were mentioned, but the focus of this paper was on development of the method. 
None of the examples mentioned, however, were composite structures. Noor and Peters (1980) 
presented a global/local approach to perform nonlinear finite element analysis applied to large de- 
flections of isotropic clamped beams and shallow arches. The discretized area was limited to a small 
number of deformation modes which were described by functions commonly used in static pertur- 
bation methods. Rayleigh- Ritz approximations were used to reduce the vectors of unknown nodal 
forces and displacements with these functions, and therefore to reduce effort required to solve the 
nonlinear analysis. 

Han and Abel (1985) developed and used an adaptive hybrid substructuring technique where the 
entire body was initially modeled with three-dimensional linear elastic superelements. As the 
loading increases and localized stresses approach a user- specified critical stress, then those areas are 
modeled with nonlinear elements. While automatic mesh refinement was implemented, provisions 
for user- interfaced refinement were made. The stiffness matrix was updated whenever a linear 
portion developed nonlinearities, and for each load step in any nonlinear portion of the model. 
Static condensation was done to eliminate degrees of freedom. While this method appears advan- 
tageous, as more of the structure develops nonlinearities its tractability may decrease. Examples 
presented were metallic structural components using the yield stress as the user-defined critical 
stress. 
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Division of a structural component into localized areas which could be considered in parallel with 
global/local analyses techniques appears to be a viable way to take advantage of increasingly more 
common and powerful parallel computers. Sun and Mao (1988) have implemented such a tech- 
nique which exploits the rapid convergence of displacement fields for finite element analyses. A 
coarse global finite element two-dimensional solution is generated, and its displacements are applied 
to the local refined models of interest. Application of the displacements is done directly or with the 
shape functions of the global model when nodal correspondence is not exact. Because quadratic 
elements were used, the local model is not constrained linearly. Each local 'zoomed' two- 
dimensional region then had a processor devoted to its solution. 
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3.0 Common Features of the Global/Local Finite 
Element Analysis Techniques 


The problems analyzed in the present study include: a) cross-ply plates under compression, b) 
angle-ply plates under compression, and c) quasi- isotropic plates under compression. The 
interlaminar stress states for a portion of the straight free edge near a central circular hole and the 
region near the hole itself were examined. Common features of the problem description and ap- 
proach for these cases and some special considerations are presented in the following section. 


3,1 Assumptions and Limitations of the Global/ Local 
Procedure 


The two-dimensional to three-dimensional global/local finite element analysis used in the present 
study has certain limitations which cannot be overcome. Some of these limitations are generally 
found within two-dimensional or three-dimensional finite element analysis, while others are attrib- 
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utable to the global/local analysis technique itself. At all times, it must be remembered that the 
finite element method is an approximate solution to a set of differential equations governing the 
elastic response of a medium. An elasticity solution, for example, would satisfy compatibility and 
equilibrium over the entire domain exactly, while the finite element method may satisfy these con- 
ditions in a piecewise manner (Reddy, 1984, 1989), (Zienkiewicz, 1986), (Kardestuncer, 1987). 
Such approximations are necessary when no exact solution is available. The limitations of the 
global/local finite element analysis are discussed later. 

The material properties used for the global and local finite element analyses are the elastic moduli 
and Poisson's ratios presented in Table 1. The U.S. Customary units values were used during the 
analyses. Most experimental procedures used to determine material properties involve testing 
laminates of a specific layup, applying a known load, and backcalculating the properties of the 
fiber/matrix system from the strain state. This approach assumes that the laminated composite is 
comprised of perfectly bonded laminae with no mechanical property variations within a particular 
lamina. These properties are smeared: the fibers and matrix are assumed to be an homogenous and 
orthotropic system. These must be converted to oriented layer-by-layer constitutive properties for 
use in determining the mechanical response of a laminated composite with a three-dimensional fi- 
nite element analysis, while classical two-dimensional plate/shell finite element analyses require a 
stiffness quantity. 


Table I. Mechanical Properties of T300/S208 Graphite/Epoxy 


E, 

e 2 

e 3 

g 12 

G,3 

g 23 

Vl2 

V|3 

V23 

19.2 Msi 

1.56 Msi 

1.56 Msi 

0.82 Msi 

0.82 Msi 

0.54 Msi 

0.239 

0.239 

0.45 

132.3 GPa 

10.7 GPa 

10.7 GPa 

5.65 GPa 

5.65 GPa 

3.72 GPa 

0.239 

0.239 

0.45 


The conversion of material properties for a system to the stiffness quantities required by the two- 
dimensional finite element analyses is done by integrating the effect of each layer through the 
thickness as done for classical lamination theory. Thus, the stiffness of different plies and their 
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orientations are combined, generally, by the constitutive processors of the finite element software 
system selected. Assumptions in classical lamination theory include: the laminate is assumed to 
be thin, made of perfectly bonded individual orthotropic laminae, each of which obeys Hooke's 
Law, and is linear elastic in its response, as laminae and as a laminate. These bonds are assumed 
to be infinitesimally thin, and not shear-deformable. The reference surface of a laminate is usually 
considered to be a laminate mid- surface, and a straight line drawn perpendicular to this surface will 
remain straight, normal and its original length throughout the deformation. All interlaminar strains 
and stresses are therefore assumed to be negligible, and all through-the-thickness variation is lost. 
Thus, the global model displacement values for a particular node are an 'average' of the actual re- 
sponse for the laminate. Because of the method to determine the the stiffness quantities used in 
two-dimensional finite element analysis, displacement or strain fields from the global model cannot, 
in general, determine effects of the individual material layers (Jones, 1975), (Reddy, 1984, 1989), 
(Zienkiewicz, 1986). 

Isoparametric elements can accurately determine the displacement field for a given converged model 
by using a polynomial description of the displacement field. Equilibrium is enforced through the 
static equations between the internal forces generated throughout the model and the applied ex- 
ternal forces. The strain field at the integration points is determined on an element-by-element basis 
as a derivative of the element displacement field. Nodal point strains, however, may be discontin- 
uous. This method of strain determination has been accepted by finite element users since the in- 
itial use of finite element analysis. However, the strains, and ultimately the stresses which are found 
by multiplying the strains by an elasticity matrix, are determined by derivatives of displacements 
and are not as accurate as when equations governing their behavior are included in the element 
stiffness matrix formulation, as in the formulation of hybrid elements. 

The technique of enforcing conditions on the strains and stresses inside an element creates a new 
class of elements called hybrid- and mixed-formulation elements, which are often problem specific, 
more expensive, and less versatile (Zienkiewicz, 1986), (Kardestuncer, 1987). With displacement- 
based elements, stresses determined along traction-free edges are not identically zero. A traction- 
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free hybrid element has been developed to overcome this difficulty (Spilker and Chou, 1980). Even 
though these elements are available, the use of displacement-based finite elements has not dimin- 
ished noticeably. If a global/local modeling technique is to be general in its application, then the 
workhorse elements used by most analysts should be used, not special elements which may be of 
limited availability. The stress field, as determined by these models, will therefore be limited in its 
accuracy and should not be taken as absolute, but rather as a strong indicator of the distribution 
and relative magnitude. 

The global/local finite element method developed in this thesis requires that displacements from the 
global model be applied to the local model in a consistent fashion. As with most computer-based 
applications, when the displacements are calculated in the global model and then applied to the 
local model, there is a limit in the accuracy involved. The global displacement field has a finite 
number of accurate digits, which are then applied to a system of equations which may encounter 
further accuracy limitations because of computer round-off errors to determine the local model 
displacement field. When. the strains are determined as gradients of these displacements, small in- 
accuracies can result in spurious strains. In addition, should repeated zooming be desired from the 
global model to the local model, and then to another local model, care must be taken that the dis- 
placements which are applied are still significant and not compromised by the displacement appli- 
cation of the previous step. The inherited errors from previous models may grow and dominate the 
solution and lead the analyst astray. 


3.2 Finite Element Software System 


All of the analyses reported herein were done on the CONVEX C220 computer in the Structural 
Mechanics Division at the NASA Langley Research Center, and all computation times reported 
are from batch jobs submitted on the CONVEX. The finite element software system used is under 
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continuing development by the Computational Mechanics Branch (CMB, formerly the Computa- 
tional Structural Mechanics Group) at NASA Langley Research Center and is called the Testbed 
(Knight, 1989), (Stewart, 1989). The Testbed is a set of multi-purpose modular processors which 
allows for the study of computational and structural analysis methods under a common framework. 
The hierarchy inherent in the Testbed is represented in the pyramid shown in Figure 1. A user 
executes or customizes procedures which interact with and invoke processors. A user of the 
Testbed requires no more detailed knowledge about the processors than would be required to use 
a commercially available software package. Installation of elements into the Testbed is relatively 
straightforward through the generic element processor (Stanley and Nour-Omid, 1990). Once im- 
plemented, benchmarking and performance evaluations of an element can be conducted. The ar- 
chitecture of the Testbed data structure is independent of the operating systems and the computer 
itself. Some processors exploit the architecture of the individual systems, but most do not. Re- 
search in structural analysis methods can take place in a cooperative environment at the procedure 
or processor level, while computational methods research often delves deeper toward the base of 
the pyramid. 

The EX97 element of processor ^ESl implemented through the generic element processor was 
used to solve all of the global models. This two-dimensional element is a nine-node assumed 
natural-coordinate strain element installed in the Testbed by Lockheed Palo Alto Research Labo- 
ratories. The local model was analyzed with classical 20-node displacement-based Serendipity ele- 
ments. A family of these elements was installed in the Testbed using the generic element processor 
by the author, and is presented in Figure 2. Mesh studies were done by Vidussoni (1988) for the 
cross-ply compression loading case with an earlier version of the Testbed, and it was determined 
that two 20-node elements, with a full 3x3x3 Gauss point rule to integrate the element stiffness 
matrices, per material layer were sufficient to capture the desired information. This was assumed 
to be correct for the other families of loading conditions considered here. 
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3.3 Global Model Generation 


The global model was initially generated using processor 'CSM1' of the Testbed (Version 1.2), 
which produces a grid of plate elements. Each plate had a total thickness of 0.1 inch and planar 
dimensions of 18x12 inches. Three different hole diameters were considered: one, three, and six 
inches. Figure 3 is a schematic of the plate with a hole and the prescribed uniform x-displacement 
of 0.1 inches at the global x=0 edge imposed as boundary conditions. This edge was restrained 
from movement in the y-direction, the z-direction, and all rotations, while the opposing edge was 
fixed for all degrees of freedom. No other displacement or rotation conditions were applied. This 
corresponds to having the loaded ends perfectly clamped as for a compression test. 

It has been shown previously that interlaminar effects are generally confined to one laminate 
thickness from a straight free edge (Pagano, 1974). A fine mesh was generated in this region with 
'CSM1.' Three elements were distributed equally by 'CSM1' across two laminate thicknesses (0.2 
in.) of width. The line two laminate thicknesses from and parallel to the straight edge along the 
maximum y-dimension of the plate thus became the global/local interface for the straight free edge 
models. A mesh study was done by Vidussoni (1988) to determine a distribution of elements in this 
zone to model accurately the interlaminar effects in cross-ply laminates under in-plane compression. 
The three rows of elements generated by 'CSM1' for all of the global models were modified to have 
the distribution used by Vidussoni (1988), as shown in Figure 4, for a portion of the six-inch hole 
global model. The row nearest the global/local boundary is one laminate thickness in width and 
provides a zone between this boundary and the area affected by the free edge of the laminate. The 
next row was two-thirds of the laminate thickness, with the remaining row being one-third of the 
laminate thickness in width. This modification was performed so that all of the models had a 
one-to-one nodal correspondence along global/local interfaces. This two-dimensional to three- 
dimensional global/local analysis method is not restricted to having such a nodal correspondence 
to succeed, but to demonstrate the validity of this modeling technique it was so limited for the 
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Compression 


along edge x =0 in. 

allow only displacement in x 
along edge x =18 in. 

allow no movement 
along edge y =0 in. 

completely free 
along edge y =12 in. 
completely free 

x =0 in., y =0 in. 

allow only displacement in x 
x =0 in., y =12 in. 

allow only displacement in x 
x =18 in., y =0 in. 

allow no movement 
x =18 in., y =12 in. 
allow no movement 

Figure 3. Boundary conditions applied to the global model for compression 
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present study. The determination of the in-plane displacement values at those nodes which do not 
exist in the global model can be done in a number of ways. Previous work has used interpolation 
splines over a critical region of the global model (Ransom, 1989), or the global element shape 
functions (Schwartz, 1981), (Kelley, 1982), (Sun and Mao, 1988). The concentration of this thesis 
is on the methodology of the extraction of accurate global two-dimensional displacement fields and 
their application in three-dimensions to the local models. 

Vidussoni (1988) also examined the interlaminar effects around holes of different sizes. For the 
plate with a six-inch hole, it was determined that if the global/local boundary was a minimum of 
five laminate thicknesses (0.5 in.) away from the hole, the local analysis could determine accurately 
the stress field around the hole. During his study the global/local interfaces remained at the same 
global nodal locations, and the smaller holes were consequently farther from the global/local inter- 
face than the larger ones. Generally a large hole in a finite width plate is considered to have a more 
severe impact on the global response than a smaller hole, and because Vidussoni (1988) had success 
with this distance from the six-inch hole, the minimum distance from the global/local interface to 
the hole edge was selected to be five laminate thicknesses for all of the hole diameters. No attempt 
was made in the present study to further minimize this distance. To provide an ample buffer zone 
between the displacement application region of the local model and the hole edge, five rings of el- 
ements of equal distribution were placed around the central circular hole for all of the hole sizes, 
as seen in Figure 5. The number of radial spokes emanating from the hole was varied to obtain 
elements with comparatively low in-plane aspect ratios and curvatures, with the larger holes having 
a higher number of spokes and a greater number of elements around the hole region. The square 
of nodes around the holes, in this figure, was identified as the global/local boundary for all lami- 
nates. This identification was done a priori, but this type of global/local analysis is not restricted 
to foreknowledge of the exact extent of the three-dimensional stress field. However, the global/local 
interface selected should not intersect the three-dimensional effect, as will be demonstrated. The 
global model was then analyzed for the angle-ply and quasi-isotropic laminates. Cross-ply lami- 
nates contain additional planes of symmetry, so that the in-plane modeling requirement is one- 
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b) 


Figure 4, 


Nodal distributions generated by 'CSM1* and as modified for the straight edge for a six-inch 
hole model: a) Original b) As modified 
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quarter of the global model. In order to maintain a one-to-one nodal correspondence along the 
global/local boundaries, additional modifications to the nodal distributions were required. The 
rings of elements around the hole region had their nodal distributions altered so four concentric 
rings of circular elements were obtained. These rings are the same as those used in the local hole 
models for all of the laminates, and are described in detail in the following section. 


3,4 Local Model Generation 


The local straight edge model had the same modified in-plane nodal distribution as the global 
model. Three 20-node elements with full integration, resulting in nine integration points in the 
y-direction, were used across the local model width. For laminates with off-axis plies, the global 
model had an in-plane nodal distribution inside the global/local boundary which was not coincident 
with the local curved edge model. To generate these distributions, the grid provided by 'CSM1' 
was modified slightly to produce concentric rings of elements with specified widths around the hole 
(Figure 5). Because interlaminar stresses are often limited to a region which is approximately equal 
to the laminate thickness, two element rings of equal width were placed within one laminate 
thickness of the hole edge. The next two rings were specified to be one and two laminate thick- 
nesses in width, respectively, along any radial line extending outward from the curved edge. The 
final ring of elements was allowed to vary in its planar dimensions in order to mate the curved el- 
ements of the hole region with the surrounding rectangular elements maintaining a minim um ele- 
ment width of 0.1 inches throughout this matching region. 

The modified in-plane nodal distribution produced by 'CSM1' was then expanded through the 
thickness for both the straight edge and the curved edge models. Half of the laminate thickness is 
sufficient for each local model due to the symmetry of the geometry and loading conditions about 
the laminate mid-plane. A total of eight elements was used through the thickness for each of the 
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Figure 5. Nodal distributions generated by *CSMl* and as modified for the curved edge region for a 
six-inch diameter hole; a) Original b) As modified 
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local models. Each local model was bounded by the prescribed global/local interfaces, appropriate 
symmetry planes, and the traction-free edges. 


3,5 Solution Procedure 


The nodal point locations, element connectivities, material property data, and stacking sequence 
are input to the Testbed. Element-level stiffness matrices are formed within the generic element 
processor and then assembled into a global stiffness matrix for a given model. With the applied 
displacements as constraints on the global matrix, a set of static, linear equations is solved. Several 
of the equation solvers available in the Testbed exploit the vector computer architecture of the 
CONVEX computer. A processor called 'ITER,' which contains several preconditioned conjugate 
gradient iterative methods, was used for each global solution, all angle-ply local models, and any 
complete three-dimensional models. Another solver called 'BAND,' which contains a direct 
Choleski solver designed for matrices stored in a novel variable-band format, was used only for the 
cross-ply local models. 'BAND' is restricted, by the amount of available incore memory for a given 
computer, in the total number of equations that can be solved and could not be used for all of the 
models. The solution rate for iterative solvers is dependent on the speed of the solver and the rate 
of convergence, both of which are dependent on the method selected, and the parameters governing 
the initial guess from which the iteration technique proceeds (Stewart, 1989), (Poole and Overman, 
1989). For the present work, no effort was made to minimize the solution time, and the first set 
of parameters for which the solution converged was used, with an incomplete Choleski precondi- 
tioned conjugate gradient (ICCG) iterative method selected. A noteworthy point is that all of these 
analyses were performed while the author was physically at Virginia Polytechnic Institute and State 
University with results transmitted, via the Internet, to university computer systems for post- 
processing. 
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3.6 Post-processing of Displacements and Stresses 


Interlaminar stresses near the straight free edge in proximity to the hole and around the hole were 
examined. All of the stress results for the straight and curved edge models were normalized by the 
absolute value of the far-field in-plane stress (<r*) as predicted by classical lamination theory for a 
plate without a hole. This value is 57,861 psi for all of the cross-ply laminates. The angle-ply 
stresses range in value due to their varying in-plane constitutive properties: 88,088 psi for 

[± 15 2n ] s , 15,821 psi for [±45 2n ] s , and 51,992 psi for [± 15 n /±45 n ]s- The quasi-isotropic normal- 
ization constant is 42,152 psi for the [0 n /±45 n /90 n ]s stacking sequence. All of the in-plane dis- 
placements are scaled by a factor of 20 for viewing, while the z-direction displacements are scaled 
by the in-plane scaling factor multiplied by ten for the local models, unless noted otherwise. 

Integration point stresses are generally considered to be most accurate (Barlow, 1976) but are often 
at inconvenient locations, thus nodal stress values were obtained by extrapolating from the inte- 
gration point stresses when presentation of the integration point values is precluded. Information 
about the extrapolation functions which were used, including their origin, derivation, and several 
sets of examples are presented in Thompson (1990). All figures which contain stress contours re- 
quired the nodal stress values, because the software system used to present this information allowed 
either elemental values from an element centroid or nodal values to be imported, not integration 
point values. Whenever multiple values for the nodal stresses occurred for those nodes connected 
to several elements, those which were known to be continuous, were input to the software, as the 
average of the values. Intralaminar stresses are continuous on an element-to-element basis in-plane 
and out-of-plane for a given material layer, but exhibit jumps in values through a laminate when 
layer properties or orientation differ. Interlaminar stresses are continuous out-of-plane, regardless 
of changes in layer characteristics through the thickness. 
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4.0 Cross-ply Laminates Under Compression 


4.1 Introduction 


The relative success of this two-dimensional to three-dimensional analysis technique will be judged 
based on its ability to determine the localized three-dimensional displacement, strain, and stress 
fields in a more efficient manner, given only information from a global two-dimensional displace- 
ment field, than a complete three-dimensional analysis. All of the laminates considered are sym- 
metric, and thus no coupling exists between the applied uniaxial compression loads and bending. 
The lamination sequences under consideration have an equal number of 0° and 90° plies, but are 
either clustered, reversed, or dispersed. Because integration through the laminate thickness of the 
material properties of the individual laminae is performed in order to deter min e the laminate 
stiffnesses required for the global finite element solution, the global displacement field for these 
stacking sequences will be indistinguishable. The influence of the different stacking sequences on 
the interlaminar response will be determined by the local analyses. The effect of stacking sequence 
change and varying hole diameter on the interlaminar stress distributions is studied. A complete 
three-dimensional finite element analysis is completed for one of the plates and used to demonstrate 
the accuracy and efficiency of the global/local analysis technique. 
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4.2 Problem Description and Approach 


A family of plates with planar dimensions of 18x12 inches and 0.1 inch thickness is under study. 
These plates consist of laminated cross-plies with a central circular hole, and the following laminate 
orientations: [02n/902n]s, [90 2 „/02„] s , and [0 n /90 n /0 n /90 n ] s with n equal to one. These laminates will 
be referred to as clustered, reversed, and dispersed, respectively. Comparison of the stress fields for 
the stacking sequences is completed for the plates with a six-inch diameter hole. In order to de- 
termine the effect of different hole diameters, three were considered: one, three, and six inches. 
For validation purposes, a complete three-dimensional analysis was completed for the [0 2n /90 2n ]s 
laminate with a six-inch diameter hole. 

Since the in-plane response of the laminates for each hole size is identical, only one global analysis 
was performed for each hole size, while the local models were considered individually. The global 
model contains four concentric circular rings of elements with an additional ring of quadrilateral 
elements within the central square region of the plate surrounding the hole, and a graded element 
distribution along the straight free edge with the maximum y-dimension. This nodal distribution 
is described in detail in the Common Features Section. Because of the in-plane loading and the 
cross-ply stacking sequences, analysis of one-quarter of the global model would have been suffi- 
cient, but the analysis was performed on the entire global model because of the way it was generated 
by the 'CSM1' processor. Table 2 contains information about the characteristics of the global 
model for each of the three hole sizes, and the undeformed grids are presented in Figure 6a), 
Figure 7a), and Figure 8a) for the plates with six-, three-, and one-inch diameter holes, respec- 
tively. Table 2 includes information about the size of the hole in the plate, the number of nodes 
in the model, the total number of degrees of freedom, the number of elements in the entire model 
and around the hole area, the number of spokes around the hole, and the total computer analysis 
time for the model in question. 


Cross-ply Laminates Under Compression 


27 



Tabic 2. Two-Dimensional Cross-Ply Compression Global Model Characteristics 


Hole 

Dia. 

(in.) 

Nodes 

Degrees of 
Freedom 

Elements 

Spokes 
Around Hole 

Global 

Solution Time 
(CPU seconds) 

Total 

Hole 

1.0 

3,208 

19,248 

768 

160 

64 

488.9 

3.0 

2,824 

16,944 

672 

200 

80 

435.2 

6.0 

2,440 

14,640 

576 

240 

96 

242.2 


A local coordinate system with its origin at the geometric center of the plate was chosen and used 
for the straight and curved edge models with an angle theta (0) defined as a rotation about the z-axis 
from the local x-axis toward the local y-axis. These axes are coincident with the global x-, y-, and 
z-symmetry planes, as can be seen in Figure 9. The laminates consist of only 0° and 90° layers; 
and therefore, two additional planes of symmetry exist through the geometric center of the plate. 
These are aligned with the local coordinate system and are identified by their local location. These 
planes affect both local areas of interest: the straight edge model extends from the local x=0 
symmetry plane to the load application edge, while the local model for the curved edge was reduced 
to one-quarter of the square region around the hole between the local x = 0 and y = 0 symmetry 
planes. It should be noted that these symmetry planes do not necessarily exist for more general 
ply orientations or loading conditions, and this must be considered before the application of this 
global/local modeling technique. A more general class of laminates is considered in the following 
chapter. 

In-plane displacement components were applied directly to each local model along the global/local 
interfaces as applied motions. Because the assumptions inherent in classical lamination theory are 
valid in this region, there was no predicted rotation or distortion of the normal lines to the mid- 
plane, and the in-plane displacements could be applied constant through the laminate half- 
thickness. The placement of the global/local interfaces in areas where the normals were not 
distorted and where the in-plane stress gradients are low was deliberate. These exact locations were 
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Figure 6. 



a) 



Global Grids for lO^/^O^nL in Compression with Six-Inch Diameter Hole: a) Undeformed 

b) Deformed 
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Figure 7. Global Grids for lO^/Wjj], in Compression with Three-Inch Diameter Hole: a) 

Undeformed b) Deformed 
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b) 


Figure 8, Global Grids for m Compression with One-Inch Diameter Hole: a) 

Undeformed b) Deformed 
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described in the Common Features Section and can be seen in Figure 9, Figure 10, and 
Figure 11. Whenever x-direction displacements (u) were applied, the x-direction rigid body trans- 
lation of the local x = 0 symmetry plane was removed. Along the local x= 0 symmetry plane for 
all of the local models, only the x-direction displacement (u=0, because of removal of the rigid 
body motion) was applied constant through the thickness. No y-direction displacements (v) were 
applied, because the local x = 0 symmetry plane passes through boundary layers and the application 
of constant y-direction displacements would inhibit the development of the interlaminar stress field. 
Along the local y = 0 symmetry plane which is a boundary for the curved edge local models, only 
the y-direction displacements from the global model were applied constant through the thickness, 
and these displacements are all zero from the global solution. No x-direction displacements were 
applied here for the same reason as no y-direction displacements were applied along the local x = 0 
symmetry plane. The local z = 0 symmetry plane is at the mid-plane of the laminate. This plane 
was fixed in the z-direction for all of the local models and was the only location throughout the 
local model with applied z-direction (w) displacements. There were no applied conditions at any 
of the nodes along any of the free edges of the local models which did not coincide with the 
global/local interfaces or the symmetry planes. The end of the full straight edge model which is 
coincident with the load application region for the global model is constrained to move only in the 
x-direction 0.05 in., after removal of the rigid body translation of the local x= 0 symmetry plane. 

In-plane loading of a cross-ply laminate results in two interlaminar stress components, a normal 
stress (< 7 r ) and a shear stress (t^), during the transition from the internal stress state to the 
traction-free edge. The existence of the interlaminar shear stress component (r^) in the three- 
dimensional stress state precludes the application of y-direction displacements through the bound- 
ary layer since no predictions can be made from the global model about the exact character of these 
displacements. Constant displacements should not be applied to any three-dimensional element in 
any direction where the displacements cannot be expected to be constant. Along straight edges in 
cross-ply laminates under in-plane compressive loading, no interlaminar shearing occurs in the x-z 
plane ( y „ = 0). Therefore, the application of constant x-direction displacements in the z-direction 
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would not violate any kinematic conditions inherent in the three-dimensional solution of the 
boundary layer. A cut was made through the boundary layer region to form the reduced local 
straight edge models with the same x-direction dimension as their corresponding local curved free 
edge models. This tremendously reduced the modeling effort for the straight free edge region. The 
displacement application procedure for these reduced local straight edge models, shown schemat- 
ically in Figure 10, is similar to the full local straight edge models except for the treatment of the 
cut itself. Since the cut is parallel to the local x = 0 symmetry plane, it has the same rule of of dis- 
placement application: only the global x-direction displacements, with the rigid body translation 
of the symmetry plane removed, are applied constant through the thickness. A comparison of the 
interlaminar stress fields between the complete three-dimensional model and the local models is 
presented in the following sections and demonstrates that using a reduced local model for this par- 
ticular loading and stacking sequence will yield accurate answers. 

The complete three-dimensional model consisted of one-quarter of the in-plane area of the plate 
and one-half of the plate thickness. The common regions modeled by the complete and local 
models were designed to have identical nodal distributions so that the stress fields could be com- 
pared at exactly the same locations. The displacement application procedure considered only the 
x-, y-. and z-symmetry planes, and the clamped plate ends (Figure 12). The undeformed three- 
dimensional meshes for the [02 n /90j n 3s laminate with the six-inch diameter hole are shown in 
Figure 13a) and b), Figure 14a), and Figure 15a). Both an unsealed and a scaled undeformed 
model are shown in Figure 13. To see the detail of the model through its thickness, its presented 
z-dimensions were multiplied by ten for the remainder of this study. The unsealed plate shows the 
relative dimensions of the model that is being analyzed: a thin laminate with a large hole. Hole- 
radius to laminate-thickness ratios are 60, 30 and 10. The length to laminate-thickness ratio is 180. 
The width to laminate-thickness ratio is 120. Information presented in Table 3 indicates the re- 
sources required for the complete three-dimensional model, and the global/local solutions for the 
full and reduced straight edge models, and the curved edge model for the [0 2n /90 2n ] s laminate with 
a six-inch diameter hole. 
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Table 3. Comparison of Complete and Global/Local Analyses for l®2n/^®2nls Laminate with Six-Inch 
Diameter Hole 


Model 

Type 

Nodes 

Degrees of 
Freedom 

Elements 

Local/Total 
Solution Time 
(CPU seconds) 

Complete 3-D 

5,759 

17,277 

1,152 

/ 

~/2,3<jfc6 

Global 2-D 

2,440 

14,640 

576 

-/ 242.2 

Full Straight 

1,405 

4,215 

240 

162.5/404.7 

Reduced Straight 

881 

2,643 

144 

99.8/342.0 

Curved 

2,559 

1 | 

7,677 

480 

436.0/67J.2 








The full straight edge local model is presented in Figure 16a). The smaller elements near the local 
x= 0 symmetry plane form the reduced local model pictured in Figure 17a). Forty percent of the 
full straight edge model contains elements away from the curved free edge of the'frole, which are 
not the focus of this study. Should the global/local analysis be successful without these elements, 
then they need not be included in future analyses. Quantitative information about -each reduced 
local model is presented in Table 4. The local x-direction was modeled with six elements for the 
local model near the six-inch diameter hole, five elements for the three-inch case, and four elements 
for the one-inch case as shown in Figure 17a), Figure 18a), and Figure 19a), respectively. Each 
of the reduced local models has a length equal to the hole radius plus five laminate thicknesses in 
the local x-direction, matching the x-direction dimension for the corresponding curved free edge 
model. Because the stress gradients away from the edge are small and increase with proximity to 
the edge, it was advantageous to grade the mesh accordingly so that the element density increased 
closer to the free edge. The grading employed was that determined by Vidussoni (1988), as de- 
scribed in the Common Features Section, modeling the local y-direction with three elements. 


Another view of the complete three-dimensional model is presented in Figure 15a). Not including 
the elements in the straight edge region, nearly half of the model is outside the area of interest 
around the hole. The undeformed local curved edge models for the three hole sizes are presented 
in Figure 20a), Figure 21a), and Figure 22a), in order of decreasing hole diameter. The relative 
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Figure 12. Complete Model Displacement Application for Cross-plies in Compression 
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a) 



b) 


Figure 13. Complete Models for lO^n/^O^nL viewing Straight Edge near Six-Inch Diameter Hole: 
a) Actual Model b) Scaled Moad 
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a) 



Figure 14. Complete Models for [02n/^®2nt v*® 1 ** 1 ^ Straight Edge near Six-Inch Diameter Hole: 
a) Undeformed b) Deformed 
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a) 



Figure 15. Complete Models for dewing Curved Edge around Six-Inch Diameter Hole: 

a) Undeformed b) Deformed 
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Table 4. Three-Dimensional Cross-Ply Compression Local Straight Edge Model Characteristics 


Hole Dia. (in.)/ 
Lay-up 

Nodes 

Degrees of 
Freedom 

Elements 

Local/Total 
Solution Time 
(CPU seconds) 

1.0/ 

[(W902n]s 

619 

1,857 

96 

70.8/559.7 

3.0/ 

[02n/902„]s 

750 

2,250 

120 

84.5/519.7 

6.0/ 

[02n/902n]s 

[902n/02n]s 

[0 n /90 n /0„/90 n ] s 

881 

2,643 

144 

99.8/342.0 

99.4/341.6 

94.2/336.4 


dimensions of the rings around each of the holes can be seen in these figures. Further information 
is provided in Table 5. In-plane curvature of the elements was kept low, with a larger number of 
spokes for the larger hole diameters. Element arc lengths along a constant radial distance from the 
hole were varied to generate radial nodal lines from the hole edge to the global/local boundary 
around the hole. Care was taken to ensure that the mid-side nodes for a given element in the re- 
vised grid were exactly at the mid-point between the comer nodes, since this is an inherent as- 
sumption behind the governing equations of the displacement-based elements. Incorrect placement 
of the nodes will affect the solution. 


4.3 Results and Discussion 


The two-dimensional global model was analyzed, and its appropriate displacements were applied 
to the three-dimensional local models. After determining the local nodal displacements, the strains 
and stresses at the integration points were calculated. Examination of solution times presented in 
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Table 5. Three-Dimensional Cross-Ply Compression Local Curved Edge Model Characteristics 


Hole Dia. (in.)/ 
Layup 

Nodes 

Degrees of 
Freedom 

Elements 

Model/Total 
Solution Time 
(CPU seconds) 

1.0 

[0 2n /90 2n ] s 

1,755 

5,265 

320 

245.4/487.6 

3.0/ 

[0jn/90j n ]s 

2,157 

6,471 

400 

326.5/568.7 

6.0/ 

[0 2n /90 2n ] s 

[90 2n /0 2n ]s 

[0„/90 n /0 n /90„]s 

2,559 

7,677 

480 

436.0/678.2 

443.0/685.2 

440.1/682.3 


Table 3 indicates that the two-dimensional to three-dimensional global/local finite element analyses 
presented were accomplished with a reasonable expenditure of computer resources. 

There is a wide range of information available in the open literature about the interlaminar stress 
state in composite structures, and an excellent review was written by Salamon (1980). Most of the 
previous investigations have concentrated on the stress state for laminates with no internal cutouts, 
under u niaxial tension or compression. Among those who have completed work for cross-ply 
laminates are (Pipes and Pagano, 1970), (Wang and Crossman, 1977). Along straight free edges 
of cross-ply laminates under in-plane loading away from the load application region, the 
interlaminar normal stress (<r t ) is maximum at the mid-plane. Considering z-locations closer to the 
top free edge, the in-plane stress character is unchanged but diminishes in magnitude. The 
interlaminar shear stress (t„) is zero at the mid-plane and top free surface, and has local maxima 
and minim a at the 0° and 90° ply interfaces. Its character is dependent on the relative positions of 
the ply orientations at each interface. The remaining interlaminar shear stress (r ra ) is zero 
throughout the laminate. At the straight free edge, the interlaminar stress state that develops for 
cross-ply laminates is dependent on the Poisson's ratio mismatch between the layers at a particular 
interface. Because the 0° and 90° plies have different Poisson's ratio values, they tend to distort in 
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Figure 21. 



b) 


Local Curved Edge Models for I°2n/ W 2nt Three-Inch Diameter Hole: a) 

Undeformed b) Deformed 
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a) 




b) 


Figure 22 . Local Curved Edge Models for |0 2n /^2i|]s ***** One-Inch Diameter Hole: a) 

Undeformed b) Deformed 
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different directions under loading. This stress component develops to enforce compatibility of the 
displacement field, while the distribution of <r* provides a countering moment so that there is no 
net moment at the edge. The 0° and 90° layers exhibit no shear coupling in the x-z plane. 

Consideration of the deformation, and hence the stress, state around the hole, however, is more 
complex. While the global cross-ply stacking sequence is fixed, subsequent angular positions 
around the hole will experience a change in the local stacking sequence. For different angular po- 
sitions around the curved edge, the 0° and 90° plies experience a different local effective stiffness. 
Thus, the local Poisson's ratio mismatch and shear coupling between the plies will change as a 
function of the angular location around the hole. This fact, coupled with the change in load around 
the hole and through the laminate half-thickness ensures that the character of the interlaminar stress 
components around holes in cross-ply laminates is complex to predict, and three-dimensional in 
nature (Dana, 1974), (Rybicki and Schmueser, 1978), (Salamon, 1980). Along the mid-plane at the 
hole edge, the interlaminar normal stress is a maximum and reduces in magnitude through the 
laminate thickness for a given angular location. The two shear stress components (t 9i , and t„) are 
both zero at the mid-plane and top of the laminate and exhibit peaks at interfaces of the 0 and 90 
layers. For the purpose of this study, only those stress components which are significant are pre- 
sented for a given interface. 


4.3.1 Two-Dimensional Global Model Results 


The deformed global grid for the plate with a six-inch diameter hole and a [0 2 „/90 2 „] s stacking se- 
quence is shown in Figure 6b), the plate with three-inch hole in Figure 7b), and the plate with 
one-inch hole in Figure 8b). All displacements were scaled by a factor of 10 for viewing. Note that 
the ends (vertical edges) of each plate remained straight and their original length as specified by the 
boundary conditions. The circular holes became noncircular, and there was a bulging of the pre- 
viously straight edges of the plates. As expected, the deformations and influence due to the hole 
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extend farther into the plate with larger holes. The line of nodes two laminate thicknesses from the 
maximum global y-dimension is the global/local interface for the straight edge models. The cur- 
vature at this location is transferred to the local model through the displacement application pro- 
cedure described earlier. The square of nodes around each hole exhibits concavity on its vertical 
sides and convexity on its horizontal sides. These are the global/local interfaces for the curved edge 
models, and its distortion will be transferred to the local curved free edge model. 


4.3.2 Three-Dimensional Local Model Results 


The straight free edge local analyses shall be presented prior to those of the curved free edge. The 
success of the global/local analysis technique shall be determined by the comparing the interlaminar 
stress fields from the full and reduced local finite element models to the complete three-dimensional 
model. Comparison of the local stress fields for different ply orientations for a plate with a single 
hole diameter and the effect of hole size for a given layup were completed. The curved free edge 
results follow the same format: verification of the global/local model results with respect to the 
complete model, the ply orientation study for the largest hole, and finally determination of the in- 
fluence of different diameter holes for one of the stacking sequences in a finite- width plate. 

4.3.2. 1 Verification of Global/ Local Modeling Technique for the Straight Edge 

If the two-dimensional to three-dimensional global/local analysis results are indistinguishable from 
those of the complete three-dimensional analysis and determined with less effort, then this tech- 
nique can be considered both verified and successful. Specific information about the different 
models is presented in Table 3 on page 37 The combined computer time to solve the global model 
and the straight and curved edge local models is less than two-thirds the complete three-dimensional 
model. However, since all of the existing planes of symmetry were not exploited for the two- 
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dimensional global model, its solution time could have been reduced further. The presented two- 
dimensional to three-dimensional global/local analysis methodology is more efficient than a 
complete three-dimensional analysis, because only those areas of interest are included in the three- 
dimensional model(s). Parallel computation schemes become attractive since the local models are 
independent of one another, but are not explored here. 

The undeformed and deformed complete models for the [0 2n /90 2n ] s laminated quarter plate with a 
six-inch diameter hole are presented in Figure 14a) and b) and Figure 15a) and b). In-plane def- 
ormations were scaled by a factor of twenty, while the z-direction displacements were multiplied 
by 200. The plate is visibly thicker because of a through the thickness Poisson effect over the entire 
plate, and most noticeably near the hole. The local full straight edge and the complete models are 
indistinguishable through visual inspection over the straight free edge region. Because of the 
boundary conditions, those elements nearest the plate ends have large deformation gradients in the 
transition from the free z-direction expansion over the whole plate to the enforced no z-direction 
expansion region at the clamped ends. Nearest the hole along the straight edge, elements are bowed 
away from the hole, as in the global model. Along the local x = 0 symmetry plane, the deformation 
state changes from z-directional expansion due to the far-field loading conditions to z-directional 
contraction because of the free edge effects. It can also be seen that this change is confined to the 
elements within one laminate thickness of the free edge, as found previously by others (Salamon, 
1980). Visual comparison between the reduced local model and the complete model is more diffi- 
cult because of the difference in scale between the models. The reduced model exhibits Poisson 
deformation and the element bowing found in the full and complete models. 

The comparison of the stress field along the straight edge as determined by the two approaches is 
encouraging. The interlaminar normal stress (<x z ) distributions for the integration points nearest the 
mid-plane and free edge for the complete model and the reduced local model are presented in 
Figure 23a) for the entire length of the plate. For this stress component, the complete model has 
extreme values at the clamped ends because the no expansion is allowed in this region. At those 
integration points common to both models, the < j z distributions are indistinguishable. The 
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interlaminar shear stress (t^) distribution at the interface between the 0° and 90° layers at the inte- 
gration points nearest the free edge for the complete three-dimensional and the reduced local model 
are shown in Figure 23b). Again, the stress values for these two models are coincident. Both of 
the models correctly determined that the interlaminar shear stress (t«) is zero for this laminate. 
The two-dimensional to three-dimensional global/local finite element analysis produced developed 
for cross-ply laminates along a straight free edge near a hole may be considered valid, verified, and 
more efficient than a complete three-dimensional finite element analysis since the same results were 
determined with less computational effort. 

The imposition of the constant x-direction and y-direction displacements from the global model 
on the local straight edge model at the global/local interfaces was an integral assumption of the 
global/local modeling technique. However, examination of the complete model displacement field 
in this area reveals that there is some variation through the half-thickness, but that this variation 
is less than 1%. The imposed global displacements on the local three-dimensional model are 
smeared through the thickness, and therefore, are not identical to the three-dimensional displace- 
ment field. However, since the global/local interface was close but outside the boundary layer, the 
local model could develop the correct three-dimensional displacement field at the free edge, and 
thus, determine the free edge stress field correctly. 

4.3. 2.2 Straight Free Edge Model Results 

The deformed reduced local six-inch hole models for the different stacking sequences are presented 
in Figure 17b), Figure 24b), and Figure 25b). Each local model is visibly compressed in the x- 
direction. Also noticeable is an in- plane bowing of the global/local interface, as expected when the 
applied displacement conditions contain such curvature effects. In Figure 17b) the elements 
modeling the 90° plies are seen to restrain those elements modeling the 0° plies from movement 
away from the center of the plate. Because of the clustering of plies with common orientation, this 
interlaminar effect is more pronounced. The free edge becomes noticeably thinner for this stacking 
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sequence. By comparing Figure 17b) and Figure 24b), the effect of reversing the order of the plies 
can be seen with the [90 2n /0 2n ] s local model having the opposite distortion at the free edge, because 
of the reversal of order of the plies at the interface. By comparing Figure 17b) and Figure 25b), 
its is determined that the dispersed plies in the [0 n /90 n /0 n /90 n ]s laminate have the same type of de- 
formation as the [0 2n /90 2n ]s plate, but reduced in magnitude, because the alternating layers temper 
the effect of the Poisson's ratio mismatch between the plies. 

From the through-the-thickness stress distribution in Figure 26a), the maximum interlaminar 
normal stresses are found at the integration points near the mid-plane of the laminate. The char- 
acteristic shape determined by others for this stress component is maintained for each x-location 
within the local model, as seen in Figure 26b). In the deformed view of the global model there is 
displacement in this area, indicating that the hole is influencing the displacement field of the previ- 
ously straight edge, and thus the stress state. This is reflected in Figure 26b) because with in- 
creasing distance from the centerline of the hole, the absolute magnitudes of the stress increase. 
For the [0 2n /90 2n ] s orientation, the values ranged from 1.4% of the far-field stress at the local x = 0 
symmetry plane to a maximum magnitude of 2.24% away from the hole. 

The through-the-thickness distribution of the interlaminar shear stress (r^) is presented in 
Figure 27a) at the z- 0.025 in. interface for each x-location in the local model. Because there are 
no integration points at the element boundaries, all interfacial stress values are an average of the 
integration point stresses on either side of the interface. From stress-reporting studies completed 
by the author, it was determined that this approach was sufficiently accurate for the nodal distrib- 
utions used here (Barlow, 1976, 1989). This stress is maximum at the interface of the 90° and 0° 
layers and zero at the mid-plane and top free edge of the laminate. As with a Zt the character of both 
its through-the-thickness and in-plane distribution becomes more exaggerated farther from the hole, 
shown in Figure 27a) and b), respectively. The range in stress values from the centerline of the 
plate to the farthest local model integration point is from 1% to 1.8% of the classical lamination 
theory in-plane normal stress, demonstrating that the effect of the large hole in proximity to the free 
edge is to reduce the absolute magnitude of the stress. Ty* has an in-plane distribution that is gen- 
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Figure 26. <r for Reduced Straight Edge Model for [02ii/^2nk : a ) Th rou fch -the -thickness b) Near 

Mid -plane 
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erally zero, but rises sharply near the free edge. As these gradients are within one quadratic element, 
further refinement or higher order elements may yield different values at the free edge. In the limit, 
however, the exact stress value is unattainable, because displacement-based elements do not satisfy 
the traction-free condition exactly (Salamon, 1980), and (Spilker and Chou, 1980). 

Figure 28a) is the through-the-thickness distribution, while Figure 28b) is the mid-plane distrib- 
ution of &i at a straight edge for the [90 2n /0 2 r»]s lamination sequence near a six-inch hole. As for 
the previous lamination sequence, the maximum values are found along the mid-plane in 
Figure 28a). The change in stacking sequence causes a change in sign of this stress component. 
However, reversing the stacking sequence does not produce a mirror image of the in-plane distrib- 
ution, seen by comparing Figure 26b) and Figure 28b). The [90 2n /0 2n ] s lay-up produces maximum 
tensile values (2.8%) which are greater in magnitude than those compressive values for the 
[0 2 n/90 2 n]s laminate (2.4%). Again, the larger absolute magnitudes are found at the farthest inte- 
gration points from the hole centerline. At the integration point nearest the centerline of the plate, 
the absolute value of this stress component is larger for the [90 2n /0 2n ]s laminate, 1.6% versus 1.4%. 
The influence of the nearby hole is to reduce the magnitude of a z along the free edge from a value 
of 2.8% to 1.6% of the in-plane normal stress in the same pattern as that for the [0 2n /90 2n ]s lami- 
nate. 

As for the interlaminar normal stress component, reversing the order of the plies in the laminate 
changes the sign of the through-the-thickness distribution of the interlaminar shear stress, as seen 
in Figure 29a). The peak values for this stress component are again at the interface between the 
0° and 90° layers. The result of reversing the lamination sequence on the in-plane distribution of 
Tyz is to produce a stress distribution that is a mirror image of that for the [0 2n /90 2n ]s sequence, with 
the same minimum and maximum absolute magnitudes as before. Once again, as the integration 
points are farther from the hole along the free edge and thus, the effect of the hole is lessened, the 
absolute magnitude of Tyj becomes larger. 
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Figure 29. t„ for Reduced Straight Edge Model for 190^/0^: a) Through-the-thickness b) At 

z*u.025 in. 
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Unlike the other cross-ply laminates, the [0 n /90 n /0 n /90 n ]s laminate has several interfaces between 
0° and 90° plies. As a result, its through -the -thickness interlaminar normal stress distribution ex- 
hibits more sign reversals and a reduction in overall magnitude than the clustered stacking sequences 
(Figure 30a). Because of the dispersal of the 0° and 90° layers, the influence of the Poisson's ratio 
mismatch between the layers is not as concentrated as with the clustered stacking sequences, and 
exists at more locations through the laminate thickness. The in-plane distribution at the integration 
points nearest the mid-plane of cr z , as shown in Figure 30b), is similar in nature to the [0 2n /90 2n ]s 
laminate, however, its magnitude is lowered. Again, the effect of the large hole in proximity to the 
free edge alters the values for the stress values but not the character of its distribution for a given 
x-location in the local model. 

The through-the-thickness interlaminar shear stress distribution, in Figure 31a), exhibits more 
variation than the clustered lamination sequences, as expected. Peaks exist at the 90°/0° (z= 0.025 
in.) interface, with valleys at the 0°/90° (z = 0.0125 in. and z = 0.0375 in.) interfaces. The fact that 
the effect of the Poisson's ratio mismatch between the different layers is dependent on the order of 
the plies about a particular interface is shown clearly in this figure. To date, the interface of interest 
in this work has been the z = 0.025 in. location, which is the only interface common to all three 
laminates. Because the order of the ply orientations from the mid-plane above the z= 0.025 in. 
location is 0° and then 90°, its in-plane distribution should resemble the [90 2n /0 2n ]s laminate and 
not the [0 2n /90 2n ] s laminate. Since the dispersal of the plies lessens the effective Poisson's ratio 
mismatch, the stress magnitudes at this location should be less than that for the clustered lamination 
sequence. Examination of Figure 31b) supports both of these statements. 

The effect of clustered, reversed, and dispersed laminate orientations along the integration points 
nearest the free edge for cross-ply laminates with a six-inch diameter hole for both of the stress 
components of interest is found in Figure 32. First the interlaminar normal stress at those inte- 
gration points nearest the mid-plane is shown in Figure 32a). Dispersal of the 0° and 90° layers in 
the [0 n /90 n /0 n /90 n ]s laminate reduces the stress values by 50% from the [0 2n /90 2n ] s laminate over 
this region. The fact that the [90 2n /0 2n ] s laminate does not have a mirror-image stress distribution 
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Figure 30. <r z for Reduced Straight Edge Model for Dispersed Lay-up: 
Near Mid-plane 


a) Through -the -thickness b) 
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of the [0 2n /90 2n ]s laminate is clearly visible in this figure, since throughout the local model near the 
free edge its magnitude is clearly larger. For all of these stacking sequences, the presence of the 
six-inch diameter hole reduces the magnitude of the stresses. Within each local model, the peak free 
edge stress is at the farthest integration point from the hole. 

The effect of altering the order of the plies on is shown in Figure 32b) for the z= 0.025 in. 
interface. As was noted in previous discussions about this component, for the [0 n /90 n /0 n /90 n ] s 
laminate will more closely resemble the [90 2n /0 2n ] s laminate than the [0 2rk /90 2n ] s laminate at this 
interface, because of the order of the plies immediately above and below it. The magnitude for the 
dispersed laminate is reduced in magnitude from its corresponding clustered laminate by nearly 
10%, which is not as large a reduction as found for the interlaminar normal stress. An exact re- 
versal in the sign of the stress distributions for the two opposing clustered sequences is shown in 
this figure. As found for the interlaminar normal stress distributions, the integration points farthest 
from the hole in the local model had the largest magnitude of stress. This suggests again that the 
six-inch diameter hole relieves the stress field in all of the laminates considered. 

While the six-inch hole has a considerable effect on the straight free edge stress distributions for all 
of the stacking sequences, the effect of different hole diameters in a finite width plate remains to 
be considered. Examination of the deformed local straight free edge meshes for the [0 2n /90 2n ]s 
laminate for the different diameter holes in Figure 17b), Figure 18b), and Figure 19b) reveals the 
trend. The influence of the hole on the local straight free edge was most pronounced near the six- 
inch diameter hole, which caused a bulging outward of the model away from the center of the plate 
(Figure 17b). This phenomenon was observed for the model of the straight edge near the three- 
inch diameter hole, but to a lesser degree (Figure 18b). As expected, there is little in-plane curva- 
ture along the edge near the plate, with a one-inch diameter hole (Figure 19b). This hole is 
one-twelfth of the plate width and is not expected to influence the previously straight free edge as 
strongly as the holes which are one-quarter (three-inch diameter), and one-half of the plate (six-inch 
di: -ter) width. The same trend was noticed earlier with the deformed views of the global models 

for the plates. 
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Figure 32, Interlaminar Stresses nearest Straight Edge for Different Laminate Orientations 
near Mid-Plane b) at z* 0.025 in. 
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The trend found for increasing hole diameter in a finite-width plate on the deformation state of the 
local models will undoubtedly be reflected in the interlaminar stress distributions, and is considered 
next. The interlaminar normal stress distributions near the mid-plane for the [0 2 n/90 2 n]s laminates 
with three- and one-inch holes is shown in Figure 33a) and b), respectively. As found for the plate 
with a six-inch diameter hole, the characteristic shape of this stress component is maintained for 
each x-location in the local model for plate, regardless of the hole diameter. As expected, the 
straight edge near the three-inch hole is affected, but not as strongly as the straight edge near the 
six-inch diameter hole. The variation in the stress, as seen in Figure 33a), ranges from 1.7% to 
nearly 2.2% of the reference stress. The model near the one-inch hole exhibits a nearly uniform 
distribution across the x-dimension of the model, and most closely resembled the distributions 
found in the literature for plates without holes (Figure 33b). The peak stress value found was 2% 
of the in-plane normal stress predicted by classical lamination theory for plates without holes. 

In Figure 34, the interlaminar shear stress distributions at the z = 0.025 in. interface are presented 
for the plates with the two smaller holes. Like the straight edge near the six-inch diameter hole, 
there is a steep rise in the stress as it nears the free edge for each of the straight edges near the smaller 
holes. From Figure 34a), it can be noted that the influence of the three-inch diameter hole on the 
straight edge alters the magnitude of the peak stress from 1.3% to 1.6% of the reference value. 
As found for the other lamination sequences and hole diameters, the largest absolute magnitude of 
stress exists along the straight edge farthest from the hole region. As seen in Figure 34b), the 
one-inch diameter hole appears to have no effect on this stress component, since its distribution is 
uniform across the x-dimension of the local model with a value of 1.5%. 

Direct comparison of o t at the integration points nearest the mid-plane and free edge for all three 
hole sizes is presented in Figure 35a) for the [0 2n /90 2n ] s laminate, while in Figure 35b), the same 
comparison is presented for at the z= 0.025 in. interface. As the diameter of the hole increases 
in a finite width plate, the effect of the hole on the stress state along the straight free edge near the 
hole increases. The presence of the nearby hole reduces the magnitudes of all of the stress com- 
ponents, both intralaminar and interlaminar, nearest the hole centerline. The largest reduction in 
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free edge stresses in proximity to the hole was found for the plate with a six-inch diameter hole, as 
was the greatest change in absolute magnitude of the stress for both of these stress components. 
Regardless of the hole diameter, the same end-shortening displacement condition was applied to 
all of the plates. The relative stiffness of each plate will vary, however, because of the reduction in 
net section and hence the ability to resist the load because of the presence of the hole. The larger 
holes will reduce the stiffness of the 18x12 inch plate to a higher degree than the one-inch diameter 
hole, and higher intralaminar and interlaminar stresses will result. Each figure reveals that the stress 
values for the local model near the larger holes overshoots that determined for the one-inch diam- 
eter hole at the end of the local models, verifying the above statements. 

4.3. 2. 3 Verification of Global j Local Modeling Technique for Curved Edge 

As with the straight edge, validation of the two-dimensional to three-dimensional global/local 
modeling technique was done by comparing the displacement and stress fields of a local model to 
those from a complete three-dimensional finite element model Information about the models is 
presented in Table 3 on page 37. For all of the characteristics presented in this table, the technique 
of analyzing a global two-dimensional model (an inefficient one, at that) and then a local three- 
dimensional model is clearly more efficient than a complete three-dimensional analysis, but the is- 
sue of accuracy remains. 

The undeformed and deformed views of the complete model of the [0 2 n/90 2n ]s laminated quarter 
plate with a six-inch diameter hole are presented in Figure 15a) and b), respectively. From this 
viewing angle, the variation in the thickness of the plate with angular position around the hole is 
readily visible. The portion of the plate near the local y = 0 symmetiy plate exhibits z-direction 
contraction of the elements near the hole edge, while near the local x = 0 symmetry plane, the ele- 
ments undergo expansion. Near the local x=0 symmetry plane, the extent of the hole effect is 
approximately 0.5 in. into the plate, which is the same distance from the boundary of the local 
curved edge model to the hole edge. In the two-dimensional global model, the compressive load 
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Figure 35. Interlaminar Stress Distributions at the Straight Free Edge for Different Holes: a) <r z 

near Mid-plane b) at Interface 
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resulted in the development of a noncircular hole, and no hint of these thickness changes existed. 
In the three-dimensional models, the in-plane change from a circular hole to a noncircular hole 
affects the out-of-plane response. The area near the hole at the local x=0 symmetry plane is 
compressed in-plane, and must bulge out-of-plane to compensate, while the area near the local y = 0 
symmetry plane expanded in-plane and contracted out-of-plane. While the in-plane loading is 
causing expansion through the laminate thickness at angular positions near 90°, the combined effect 
of the local strain and Poisson's ratio mismatch is to reduce this growth, as seen in Figure 15b), 
along the both local axes for the two rows of elements within one laminate thickness of the edge. 
The same influence is seen for the curved edge local model in Figure 20b). Because of the differ- 
ence in the scaling between these two figures, these appear erroneously to be more pronounced for 
the local model. 

Comparison of the stress fields for the curved edge model and the complete three-dimensional 
model was completed for the interlaminar normal stress near the laminate mid-plane, and the 
interlaminar shear stress components (r 9z and r n ) at the z= 0.025 in. interface of the laminate, and 
is presented in Figure 36 and Figure 37a) and b), respectively. The stresses are identical in value 
for these models for the interlaminar normal stresses outside the region between 0 = 35° and 
0 = 50°. The difference in the stresses in this region is small, however, less than 0.25% of the in- 
plane normal stress as determined by classical lamination theory, but its presence casts a shadow 
on the global/local analysis method. Examination of the undeformed two- and three-dimensional 
models showed that between the global/local boundary and the rings of circular elements around 
the hole, there are elements which are long, tapered, and skewed. While displacement-based ele- 
ments generally determine displacement fields accurately, whenever distorted elements are used, the 
displacement field can degrade rapidly (Reddy, 1989). Thus, the two-dimensional displacements 
from 'questionable' elements in the global model, are applied to 'questionable' three-dimensional 
local elements. The complete three-dimensional finite element model contains these element ge- 
ometries as well, but in contrast to the local model, well-behaved elements surround these elements 
so that their collective stiffness buffers the effect of the distorted elements. Examination of the in- 


Cross-ply Laminates Under Compression 


74 


plane displacements in this region between those applied to the local model and those determined 
by the three-dimensional model, revealed slight variations along the global/local interfaces. The 
greatest variation was found near the same angular locations where the stress states for the local and 
complete models did not coincide, precluding the determination of the identical stress field in this 
area. The fact that local model was able to determine the three-dimensional stress field identically 
in the other regions around the hole is encouraging; there appears a certain robustness in the 
global/local analysis technique, but the application of displacements that are not necessarily accu- 
rate in the global sense will adversely affect the local solution. 

Examination of the interlaminar shear stresses for the local and complete models at the z = 0.025 
in. interface is more encouraging. The r 6z distribution for both of these models is presented in 
Figure 37a) and is coincident. The global/local analysis technique determined the character and 
magnitude of this stress component correctly for all of the angular positions. Note that magnitude 
of this component is large in comparison to the normalization stress. This is because of the 
changing local lamination sequence, and hence local material properties in the laminae, and the load 
path around the hole. The remaining stress component (? n ), in Figure 37b), is smaller in magni- 
tude than the other stress components. While this stress component reaches a value of zero at the 
hole edge, it is presented here at the same radial distance from the hole as the nearest integration 
points, which is in an area of large stress gradients. A slight difference between the complete and 
local three-dimensional model results was found in the same region as found for the interlaminar 
normal stress, which was less than 0. 1 % of the normalization stress. It is acknowledged that a wiser 
element selection process could have improved the global/local results, but the relative success of 
the global/local procedure is considered to be excellent. Within limits, the local model can com- 
pensate for the fact that the displacements are accurate, or nearly so, in the global sense but not 
locally accurate. The two-dimensional to three-dimensional global/local procedure used to deter- 
mine the stress state around a curved edge for cross-ply laminate in compression can be considered 
a valid and viable modeling technique. 
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Figure 37. Interlaminar Streaaea for Complete and Reduced Local Model at Curved Edee: a) t at 

z- 0.025 in. b) at z« 0.025 in. #2 
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4.3.2.4 Curved Free Edge Model Results 


The undeformed and deformed local curved edge meshes for the different stacking sequences around 
a six-inch diameter hole are shown in Figure 20a) and b), Figure 38a) and b), and Figure 39a) and 
b). Each of the local models are compressed in the x-direction and the holes became noncircular, 
as it did for the global models. There is a noticeable nonuniform thickening of the local models 
around the hole. Each local model, regardless of the stacking sequence, exhibited thickening near 
the local x = 0 symmetry plane and contraction near the local y = 0 symmetry plane, which is at- 
tributable to the applied load on the plate. The effect of the laminate stacking sequence manifested 
itself in the rows of elements nearest the curved free edge. For the [(hn/902 n ]s laminate, those rows 
contracted more toward the mid-plane of the laminate than the [0„/90 n /0 n /90 n ]s laminate near the 
local x= 0 symmetry plane. Because the [0 n /90„/0 n /90 n ] s plate has an alternating ply sequence, 
some of the rows of elements modeling the plies contracted while others expanded, and the effect 
of the local Poisson's ratio mismatch and shear coupling was not as concentrated (Figure 20b) and 
Figure 39b). In contrast, the [90 2n /0 2n ] s model had regions of noticeable element expansion near 
the local x=0 and local y=0 symmetry planes away from the laminate mid-plane (Figure 38). 
From the viewing angle selected, it is difficult to see that there is a difference in the hole diameter 
through the laminate half-thickness for each of these stacking sequences, but the existence of r n 
inside the hole edge suggests that hole diameter reduction is not constant. It is obvious that the 
deformation states between the different lamination sequences vary and not necessarily in an intui- 
tive manner. 

The displacements were post-processed to determine the stress fields around the hole near the 
mid-plane and at the interface of the plies. For the integration points nearest the curved edge, 
through-the-thickness stress distributions were generated for all lamination sequences and 
interlaminar stress components and are considered in detail. No distributions parallel to the mid- 
plane of the laminate or those in the r-z plane for a given angular position are presented in the in- 
terest of brevity, not as a limitation of the global/local modeling technique. 
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Figure 38. Local Curved Edge Models for 190^/0^ around Six-Inch Diameter Hole: a) 

Undeformed b) Deformed * 
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The interlaminar normal stress distribution around the hole and through the laminate thickness is 
presented in Figure 40a), while the stress contours superimposed on a deformed view of the local 
model are in Figure 41. The stress contours were produced from stress values extrapolated to the 
nodes from the integration point stresses by the extrapolation technique described in detail in 
Thompson (1990). From these figures, it can be determined that a t is largest and negative along 
the mid-plane and diminishes in magnitude through the laminate, becoming zero at the top free 
edge. At the interface between the 90° and 0° plies, the stress gradients change rather abruptly, and 
diminish more rapidly. At approximately 0 = 72°, its maximum is attained, and there are consid- 
erable gradients both through the thickness and in the angular direction. The existence of high in- 
plane displacement gradients in this region coupled with the the local Poisson's ratio mismatch 
between the layers affects the character and the magnitude of the stress component. In Figure 41, 
the stress contours are seen to be concentrated in those angular positions nearest the local x = 0 
symmetry plane and within one laminate thickness of the curved free edge. The negative values of 
a z are found in the region of the curved edge model that is contracted toward the mid-plane, and 
are concentrated in the 90° ply. Those elements in the 0° ply, while contracted toward the mid- 
plane, experience no interlaminar normal stress and appear to follow those elements within the 90° 
ply toward the mid-plane. 

The interlaminar shear stress (t»,) through-the-thickness stress distribution and contours are found 
in Figure 40b), and Figure 42. With Figure 40b), it was determined that consideration of this 
component is significant only at the interface between the plies. At 0 — 75°, this shear stress reaches 
a maximum and at 8= 10° a minimum. The growth and decay of this component through the 
laminate is immediate, even, and centered about a particular angular location. This is due to the 
larger gradients in the in-plane stress in this region caused by the redistribution of the deformation 
and stress state around the hole: greater interlaminar stresses must exist to enforce the displacement 
boundary conditions. This stress component indicates the degree of slippage of the elements past 
one another around the hole. Those elements with large positive stress values moving in one di- 
rection and those with negative values in opposition. Examination of the stress contours on the 
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Figure 40. Interlaminar Stress Distributions Curved Edge Model for * round Six-Inch Di- 
ameter Hole: Through -the -Thickness a) a z b) t #i c) 
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Figure 41. <r z Contours for Curved Edge Model for around Six-Inch Diameter Hole 
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deformed model does not reveal much information with the selected viewing angle, however, the 
elements in the 90° ply are noticeably moving upward and away from the mid-plane, while those 
in the 0° ply are slipping toward the mid-plane in the region between 0 = 60° and 8 = 90°. 

Through-the-thickness distributions and stress contours around the hole for the the remaining 
interlaminar stress component (t„) are in Figure 40c), and Figure 43. Its distribution exhibits large 
gradients near the interface. Unlike the other shear stress component, large opposing gradients exist 
on either side of the 0°/90° interface. In Figure 40c), the concentration of this stress nearer to the 
local x = 0 symmetry plane is noticed. Again, the region of the hole with a larger normal stress and 
greater degree of shear coupling is where the stresses are larger in magnitude. In Figure 43 along 
the local x = 0 and y = 0 symmetry planes, these stresses are shown to be nonzero only within one 
laminate thickness of the curved free edge. The 90° and 0° plies will tend to distort the hole in a 
different direction under the same load, and since those elements at the ply interface are constrained 
to remain in contact, interlaminar stresses arise to enforce this and generate gradients on either side 
of the interface. 

Reversing the stacking sequence of the laminate has a pronounced effect on the interlaminar re- 
sponse around the hole. Presented in Figure 44a) is the through-the-thickness distribution of the 
interlaminar normal stress at the integration points nearest the hole for a [90 2n /0 2 „] s plate with a 
six-inch diameter hole. Comparing this distribution with that for the [0 2n /90 2n ]s laminate in 
Figure 40a), it was found that these distributions are not mirror images as an inexperienced analyst 
might assume. The decay of the magnitude of <r z through the laminate is expected. At the locations 
nearing the local x = 0 symmetry plane, the mid-plane magnitude of a t for the [0 2n /90 2n ] s laminate 
is considerably reduced from its peak at 9 = 72°, while for the [90 2n /0 2n ] s laminate there is no sig- 
nificant reduction in value. Those elements nearest the hole edge for this laminate would be ex- 
pected to move in the opposite fashion to those for the [0 2n /90 2n ]s laminate, and comparison of 
Figure 20b) and Figure 38b) shows this to be true. For all of the angular positions around the 
hole, the interlaminar stress (<*,) is negative for the [0 2n /90 2n ] s laminate and positive for the 
[90 2n /0 2n ] s laminate. Those elements within one laminate thickness of the curved free edge for the 
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[ 0 2 n/ 902 n]s laminate contracted toward the laminate, in the opposite direction to those elements 
under the influence of the global loading alone. Those same elements for the [90 2n /0 2n ]s laminate 
expanded away from the mid-plane, further increasing the thickness of the plate. 

In Figure 44b), the through-the-thickness distribution of r* z is presented for [90 2n /0 2n ]s stacking 
sequence. As for the [0 2n /90 2n ] s laminate, this stress component changes shape and value smoothly 
through the laminate thickness and with angular position for a [90 2 n/0 2n ]s laminate, and, unlike the 
other stress component, is opposite in sign. Scrutinizing Figure 20b) and Figure 38b) closely 
reveals that the element boundaries immediately above and below the interface for these two 
models, which were initially perpendicular to the interface, have opposite slopes in their deformed 
state. The layers at the interface are sliding past each other in these two laminates, but in the op- 
posite direction. 

The interlaminar shear stress (r n ), is complex in nature. Large opposing gradients on either side 
of the interface are noticeable for both laminates (Figure 44c). From the difference in the sign of 
this stress for the stacking sequence reversal, the hole diameter changes would be expected to be 
different through the laminate thickness. 

When the plies are dispersed, the interlaminar stresses were found to be lower in magnitude with a 
more complex distribution through the laminate thickness than for the straight free edge models. 
The curved edge models reflect this fact. All of the interlaminar stress distributions through the 
laminate thickness are presented in Figure 45. The interlaminar normal stress, in Figure 45a), for 
this laminate exhibits negatives values both at the mid-plane and at the interface at z= 0.025 in. 
As for the [0 2n /90 2n ]s laminate, the [0 n /90 n /0 n /90 n ] s laminate has its maximum magnitude near 
8 = 72°. The 0° layers have low (near zero) a z values through the thickness of the laminate. The 
character of o z for [0 n /90 n /0 n /90 n ] s laminate around the hole at the mid-plane is similar to the 
[0 2n /90 2n ] s laminate, but reduced in magnitude. 
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Figure 45. Stresses for Curved Edge Model for Dispersed Lay-up with Six-Inch Diameter Hole: 

Through -the -Thickness a) <r 2 b) t #z c) r n 
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Due to the dispersal of the 0° and 90° layers, two distinctive peaks can be found at z = 0.0125 in. 
and z = 0.0375 in. for the interlaminar shear stress (t 9z ) with an obscured valley at the z = 0.025 in. 
interface (Figure 45b). Those 0°/90° interfaces exhibit peaks similar to the [0 2n /90 2n ] s laminate, 
while the one with the 9070° interface, as expected, has valleys like the [90 2 „/0 2n ] s laminate. Re- 
duction in the magnitudes of this stress component are noticeable for each of the interfaces, in 
comparison to the clustered laminate, while the deformed view of the local model in Figure 39b) 
reveals that the element boundaries, which were originally perpendicular to the interfaces, have 
different slopes at the interfaces. 

The remaining interlaminar shear stress (x rz ) has a more complex through-the-thickness stress dis- 
tribution around the hole because of the larger number of interfaces through the thickness than the 
other laminates studied. The peak positive value is lowered slightly for this laminate, when com- 
pared to that for the [90 2n /0 2n ] s laminate, while there is little difference in the peak negative value 
from Figure 45c). Hole diameter change through the laminate thickness for [0„/90 n /0 n /90„] s is 
lower in magnitude but with several sign reversals. 

Direct comparison of the interlaminar normal stress for changing stacking sequence at those inte- 
gration points near the mid-plane is presented in Figure 46. While the stress distributions for the 
[0 2n /90j n ]s and [ 902 n/ 0 2n ]s laminates are opposite in sign, their absolute magnitudes are not identical. 
There is a reduced peak absolute magnitude for [90 2n /0 2n ] s at a slightly larger value of 6 than for the 
[0 2n /90 3n ]s laminate. In addition, while the stresses nearest the local x= 0 symmetry plate decay in 
magnitude from the peak stress for the [0 2n /90 2n ] s laminate, this laminate does not have the same 
degree of decay. As expected, the dispersed laminate stress distribution is reduced in magnitude, 
but with the same character as the [0 2n /90 2 „] s laminate for this stress component near the laminate 
mid-plane. 

Since the only interface common to all three la min ates is at z= 0.025 in., direction comparison 
between the laminates is determined at that location for the interlaminar shear stresses (r tf2 and 
r„). These are presented in Figure 47a) and b). For t 8z , reversing the stacking sequence produces 
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an exact reversal in sign of the stresses. At z = 0.025 in. the 90° ply is above the interface, while the 
0° ply is below for the dispersed laminate, so the stress distribution is expected to resemble the 
[90 2n /0 2n ]s laminate, but with reduced magnitudes. This reduction is found to be nearly two-thirds 
at the peak stress location as presented in Figure 47a). 

Figure 47b) contains the information for the three lamination sequences about the remaining 
interlaminar shear stress (t„) at the z = 0.025 in. interface. Recall that these stresses were originally 
calculated at the integration points on either side of the interface, and that average of these values 
are presented. As for the other interlaminar shear stress, the reversal in the stacking sequence 
produces a mirror image distribution about T n — 0. While the dispersal of plies reduced the mag- 
nitude of the distribution of the [90 2 n/0 2n ]s laminate, unlike the other shear stress, its character was 
not identical. for the [0 n /90 n /0 n /90 n ] s laminate continued to increase in magnitude at the larger 
angular positions until it reached [90 2n /0 2n ]s laminate distribution. 

For the different hole sizes, presentation of the through-the-thickness distributions of the stress 
components will be omitted since they closely resemble those for the six-inch diameter hole. In- 
stead, direct comparison of the stresses is presented next. The interlaminar normal stress near the 
mid-plane of the laminate is presented in Figure 48 for the [02 n /90 2n ]s stacking sequence and all 
three hole diameters. The location of the largest magnitude stress moves toward the local x = 0 
symmetry plane with increasing hole size. The six-inch hole exhibits a lower stress magnitude than 
either of the hole cases, with the one-inch hole largest. Stress relief due to hole size in a finite width 
plate is believed to cause this phenomenon. Since the straight free edge does not strongly resist the 
outward bowing because of the hole, the in-plane stress components around the hole are redistrib- 
uted, and hence, lower interlaminar stresses result. While there is some stress relief with the three- 
inch hole, there is none for the one-inch hole, which explains the shift in peak stress location and 
value. The lack of stress relief for the one-inch hole appears to further manifest itself in this dis- 
tribution since there are both positive and negative stresses, while those for the other holes are all 
negative. The one-inch diameter hole is unable to resist the out-of-plane expansion of the plies at 
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Figure 46. Stresses for Curved Edge Model for Different Stacking Sequences Around Six-Inch 
Hole: a) <r z at Mid-plane 
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the local x = 0 symmetry plane through the interlaminar strain state. A closer examination of the 
deformed model for this case corroborates this. 

The influence of decreasing hole diameter in a finite width plate on the interlaminar shear stresses 
is presented at at the z = 0.025 in. interface in Figure 49. From Figure 49a), all three hole diam- 
eters exhibit the same character for t 9z , with the location of the maximum absolute value of the 
stress moving closer to the 6 = 90° as the hole diameter increases. The maximum magnitude for the 
three-inch diameter hole model is slightly greater than the other holes, and this suggests that there 
is a critical hole size in the laminate in a finite width plate. The results for the remaining 
interlaminar stress component, (t, z ) are found in Figure 49b). The peak locations move toward the 
local x = 0 symmetry plane with increasing hole diameter. The three-inch diameter hole is larger in 
absolute magnitude than the six-inch diameter hole for all angular positions except those greater 
than 75° and less than 90°. Once again the one-inch diameter hole exhibits both positive and 
negative shearing values suggesting that there is no stress relief for this hole size. Indeed, its char- 
acter does not resemble the larger holes as closely as the other stresses. Between 9 = 20° and 70 , 
its absolute magnitude is larger than other holes. 


4.4 Conclusions for In- Plane Compression of Cross- Ply 
Laminates 

A number of conclusions can be drawn from the analyses discussed, both relative to the state of 
stress around large holes in thin laminates and to the application of two-dimensional to three- 
dimensional global/local stress analysis in laminated cross-ply composites. The global/local con- 
clusions will be presented in a later following the development and application of a global/local 
technique on laminates with off-axis plies. 
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Figure 48. for l®2a/^®2alf nw Mid-plane Around Different Diameter Holes 
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Figure 49. Interlaminar Shear for 10^/90^1. at 2 " 0.025 in. Interface for Different Diameter 
Holes: a) %2 b) 
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Interlaminar effects at straight free edges for plates without holes have been considered extensively 
in the past. A comparative study done by Vidussoni (1988) on the effect of hole size on the 
interlaminar stress state at the straight free edge and around a hole for plates with holes was ham- 
pered by a lack of computer resources. Both the effect of stacking sequence for a large hole at the 
free edge and hole size for a particular laminate were studied in this thesis in detail. 

When considering different plates with large holes, it was found that dispersing the 0° and 90° layers 
in cross-ply laminates lowered the overall magnitude of the interlaminar normal and shear stresses, 
but did not alter the overall shape of the stress distribution for a particular interface along the 
straight free edge. At the mid-plane the [0 n /90 n /0n/90 n ]s laminate resembled the [0 2n /90 2 n]s lami- 
nate, as expected, while the [90 2n /02n]s distribution was not an exact reversal in sign for the 
interlaminar normal stress. At the z= 0.025 in. interface, the staggered laminate has a 90°/0° 
interface and thus its shear stress distribution at this interface was like the [902 n /0;>n]s laminate and 
reduced in magnitude. Reversing the stacking sequence produced an exact reversal in sign of the 
interlaminar shear stress at the interface. 

At the straight free edge, the effect of increasing hole diameter reduced the magnitude of thv stress, 
but not its in-plane distribution for a particular z-location. The six-inch diameter hole in a 12-in. 
wide plate exhibits a large stress reduction at the centerline of the hole, but increases the 
interlaminar stresses at the far edge of the local model. The one-inch diameter hole exhibited no 
effect on the straight edge nearby, while the influence attributable to the three -inch hole was be- 
tween those of the other hole diameters. 

The investigation of the curved edge region was also completed by Vidussoni (1988) for both 
stacking sequence and hole size, but concern was expressed about the computer resources. Con- 
siderable refinement around the hole was done for the present study, and the global/local interface 
was moved closer to the curved free edge for plates with smaller holes. A stacking sequence study 
for the largest hole size, and a hole size study for one of the laminates was completed. Thickening 
of laminates with holes in compression was not constant around the hole, for all of the stacking 


Cross-ply Laminates Under Compression 


97 



sequences. The degree of thickening varied with layup, hole diameter, and position around the 
hole. For applications when laminate thickness is constrained, as in the case of rivets which are set 
into the laminate, this is a source of potential damage. 

Reversal of stacking sequence for interlaminar stress components produced a reversal in sign over 
the laminate, but whether this is an exact reversal is dependent on the location in the laminate and 
the stress component in question. The interlaminar normal stress at the integration points nearest 
the mid-plane were not exactly reversed. The interlaminar shear stresses were a mirror image for 
the z = 0.025 in. interface. While "dispersion' of plies with different ply angles reduces the absolute 
value of certain stress components, its distribution is directly dependent on the order of the layers 
about the z-location in question. 

For the interlaminar normal stress around holes near the mid-plane, the effect of hole size is clear. 
Stress values around the hole for large holes are reduced in magnitude, and peak stresses are shifted 
in location. Stress relief at the straight free edge is believed to be responsible, because the one-inch 
diameter hole which has no effect on the stress field at the straight edge, experienced both a positive 
and negative stress around the hole. There appears to be a "critical” hole size in this laminate, where 
the interlaminar shear stress (t» z ) around the hole is maximized, as evidenced by the larger stresses 
for a three-inch diameter hole as compared to a one- or a six-inch diameter hole in the laminates 
presented. The hole effect for the remaining interlaminar shear stress (r ra ) is less clear. Examination 
of three-dimensional results indicates that the reduction in hole diameter for laminates loaded in 
compression is not constant for all plies. Should the hole contain a fastener, this could result in 
bearing failures and fiber crushing in the plies tending to have larger displacements. 
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5.0 Angle-ply and Quasi-isotropic Laminates 
Under Compression 


5,1 Introduction 


The global/local analysis technique discussed and presented for symmetric cross-ply laminates in 
the previous chapter is a specific implementation of a general analysis tool. Development and 
demonstration of this method as applied to laminates with off-axis plies, e.g. angle-ply and quasi- 
isotropic laminates, will be presented in this chapter. The lamination sequences considered include 
several [±©] s laminates, one [iQj/iQjjs laminate and one quasi-isotropic laminate. Several of 
the assumptions used in the global/local procedure for the cross-ply laminates are invalid for off-axis 
ply laminates. There are fewer symmetry planes to exploit, and the global strain state affects the 
local model geometry differently than the cross-ply laminates. An example of erroneously modeling 
the straight edge is compared to a correct local straight edge model to demonstrate the need to 
understand this global/local two-dimensional to three-dimensional finite element analysis procedure 
and its limitations. The interlaminar response for the different stacking sequences is then deter- 
mined for a finite-width plate with the large hole, followed by a study of the influence of the size 
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of a hole in a finite-width plate. Attempts to verify this technique with a complete three- 
dimensional analysis were unsuccessful, because the size of the resulting model was intractable. 
Verification of the technique is done on a qualitative basis by comparing the interlaminar response 
determined with this technique to that determined by other methods for laminates with off-axis 
plies. 


5.2 Problem Description and Approach 


The plates under consideration include laminated angle-ply and quasi-isotropic laminates with a 
central circular hole, and the following ply orientations: [db 15 2n ] s , [±45 2n ]s, [± 15 n /:t45 n ]s, and 
[0 n /±45 n /90 n ]s with n equal to one. All were 18 in. long, 12 in. wide and 0.1 in. thick. An inves- 
tigation of the interlaminar stress fields for a constant width plate with different hole diameters was 
performed for the [± 15 2n ]s laminate. Three different hole diameters were considered: six, three, 
and one inches. The [ ± 1 5 2n ]$ laminate with a six-inch diameter hole was modeled with a complete 
three-dimensional model and also used to demonstrate the correct and incorrect way to apply this 
global/local analysis method to the straight free edge region. 

Each plate was analyzed with a global model which contains five equal-sized rings of elements 
within the central square region of the plate surrounding the hole, and a graded element distribution 
along the straight edge with the maximum y-dimension. This has been described in detail in the 
Common Features Section. Table 6 contains information about the characteristics of each global 
model, including the diameter of the hole, the number of nodes, the total number of degrees of 
freedom, the number of elements in the entire model and around the hole area, the number of radial 
lines of nodes around the hole (spokes), and the total computational time for the model in question. 
The undeformed and deformed global grids for the [=b 15 2n ]s stacking sequence are presented in 
Figure 50a) and b), Figure 51a) and b), and Figure 52a) and b) for the plates with six-, three-, and 
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one-inch diameter holes, respectively. The undeformed and deformed global grids for the 
[±45 2n ]s, [± 15 n /±45 n ]s, and the [0 n /± 45 n /90 n ] s laminates with a six-inch diameter hole are in 
Figure 53a) and b), Figure 54a) and b), and Figure 55a) and b), respectively. 

Table 6. Two-Dimensional Angle-Ply and Quasi-isotropic Compression Global Model Characteristics 


Hole Dia. (in.)/ 

Nodes 

Degrees of 

Elements 

Spokes 

Global 

Layup 


Freedom 

Total 

Hole 

Around Hole 

Solution Time 
(CPU Seconds) 

1.0/ 

[±15 2n ]s 

3,208 

19,248 

768 

160 

64 

439.9 

3.0/ 

[±15 2n ] s 

2,824 

16,944 

672 

200 

80 

437.7 

6.0/ 

[±15 2n ] s 

2,440 

14,640 

576 

240 

96 

304.0 

[±45 2n ] s 






351.0 

[± 15„/±45 n ]s 






309.3 

[0 n /±45 n /90„] s 






429.2 


A local coordinate system with its origin at the center of the plate was chosen and used for all the 
two-dimensional and three-dimensional models with an angle (0) defined as a rotation about the 
z-axis from the local x-axis toward the y-axis (Figure 56). The local z = 0 plane was coincident 
with the laminate mid-plane. Because of the symmetry of the lamination sequences and the applied 
global loading about the mid-plane, this is a plane of symmetry for all models. The full straight 
edge model extends between both ends of the plate, while the elements bounded by the rectangular 
elements surrounding the hole form the local curved edge model, as shown in Figure 56 and 
Figure 57. 

The displacement application procedures for these local models are different from the cross-ply 
models. There are no local x = 0 and local y = 0 symmetry planes, so the ends of the local straight 
edge model are coincident with the load application regions for the global model for laminates with 
off-axis plies. These ends are clamped and undergo compressive loading by displacing the global 
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a) 



b) 


Figure 50. Global Grids for [± l5 2n L in Compression with Six-inch Diameter Hole: a) Undeformed 

b) Deformed 
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b) 


Figure 51. Global Grids for [±15^^ in Compression with Three-Inch Diameter Hole: a) 

Undeformed b) Deformed* 
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x = 0 in. end toward the global x= 18 in. end by 0.01 in. as in the global model. Along the 
global/local interfaces, the in-plane global displacement components were applied constant through 
the thickness directly to each local straight edge model. Like the cross-ply local models, the in- 
plane stress gradients were low and there was no expected rotation or distortion of the normal lines 
to the mid-plane in the regions where the global/local interfaces were located. These locations were 
described in the Common Features Section and are identified in Figure 56 and Figure 57 for both 
areas of interest. The local z = 0 symmetry plane was fixed in the z-direction for all of the models. 
There were no applied conditions at any nodes along any of the free edges of the local models which 
did not coincide with the global/local boundaries or the local z = 0 symmetry plane. 

For cross-ply laminates under in-plane loading, the only interlaminar stress components which arise 
during the transition from the internal stress state to the traction-free straight edge are the 
interlaminar normal stress (o z ) and interlaminar shear stress (t^). With no shear in the x-z plane, 
(y«= 0), the application of constant x-displacements through the thickness in the boundary layer 
region would not violate any three-dimensional kinematic conditions. Therefore, an arbitrarily lo- 
cated cut through the boundary layer could be made. This allowed the analyst to limit the local 
models to only the region of interest and is not true for any laminates with off-axis plies, like 
angle-ply or quasi-isotropic laminates, or any loading situation which would cause shearing in a 
plane where constant displacements are applied. The error in arbitrarily locating a global/local 
interface through such a region to create a reduced local model is demonstrated for a [ ± 1 5 2n ]s 
laminate with a six-inch diameter hole. 

As with the cross-ply models, the reduced straight edge model was designed to have the same x- 
direction dimension as its corresponding local curved edge model. The displacement application 
procedure for these reduced models is similar to the full straight edge models, except for the treat- 
ment of the interfaces cutting through the boundary layer. The procedure is presented schematically 
in Figure 58. This is the same displacement application procedure used for the cuts through the 
boundary layer for the cross-ply laminates in compression. Since the two-dimensional global model 
yields no information about the through-the-thickness behavior of the laminate, only an averaged 
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response at the laminate mid-plane and to enforce the clamped end condition, these displacements 
must be applied constant through the thickness. At each end of the reduced local model, planes 
of nodes pass through the boundary layer and the applied x-direction displacements were constant 
through the thickness. On the other hand, y-direction displacements were not applied along the 
cut, to avoid inhibiting the three-dimensional interlaminar effect. A comparison of the displace- 
ments for the full (correct) and reduced (incorrect) local straight edge models is presented in the 
Results and Discussion section. 

Because of the existence of off-axis plies, the entire in-plane area of the plate must be modeled 
globally, since the only plane of symmetry is at the laminate mid-plane. This was also true for the 
complete three-dimensioned model. To allow stress field comparison at the same locations, the 
common region modeled by the complete and local models had identical nodal distributions. Fig- 
ure 59 is a schematic representation of the complete three-dimensional model displacement appli- 
cation procedure. Presented in Table 7 is the computer resources required for the complete 
three-dimensional model, and the global/local solutions for the full straight edge, the incorrect 
straight edge, and the curved edge for the [± 152n]s laminate with a six-inch diameter hole. 
Table 7 contains the information that quantitatively demonstrates that the global/local modeling 
technique will undoubtedly be more efficient than a complete three-dimensional analysis, since 
there are nearly two-thirds more nodes in it than the largest local model. Its solution is currently 
intractable, and the global/local modeling procedure provides an attractive alternative. 
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Table 7. Comparison of Complete and Global/Local Analyses for |± Laminate with a Six-Inch 

Diameter Hole 



Nodes 

Degrees of 
Freedom 

Elements 

Local/Total 
Solution Time 
(CPU seconds) 

Complete 3-D 

27,112 

81,336 

4,608 

--/ + 

Global 2-D 

2,440 

14,640 

576 

304.0/- 

Reduced Straight 

1,667 

5,001 

288 

223.0/527.0 

Straight 

2,715 

8,145 

480 

344.6/648.6 

Curved 

9,648 

28,944 

1,920 

2,330.3/2,634.3 


+ model was not solved satisfactorily 


The reduced and full straight edge undeformed and deformed models are shown for the [ ± 1 5 2n ]s 
laminate with a six-inch diameter hole in Figure 60a) and b), and Figure 61a) and b). The full 
straight edge model is larger, with 65% more nodes than the reduced straight edge model. 
Figure 62a) and b), and Figure 63a) and b) are the undeformed and deformed local straight edge 
models for the [± 15 2n ]s laminates with three- and one-inch diameter holes, respectively. Since all 
of the plates had the same global dimensions, these full local straight edge models had the same x-, 
y- and z-dimensions, while the number of nodes changed as a function of the diameter of the hole. 
The plates with the larger holes had a had a greater number of nodes around the hole, which af- 
fected the nodal distribution along the straight edge in proximity to the hole. Visual inspection of 
the local straight edge models reveals that greater number of nodes in the models near the smaller 
holes, since more nodes were needed between the finely discretized hole region and the ends of the 
plate. This information is presented quantitatively in Table 8. 
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Free surface (straight edge) 


Free surface (hole) 


Figure 59. Complete Model Displacement Application for Laminates with OfT-Axis Plies in Com- 
pression 
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Figure 60. Reduced Straight Edge Models for I^IS^]^ near Six-Inch Diameter Hole: a) 

Undeformed b) Deformed 
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Table 8. Three-Dimensional Angle-Ply and Quasi-isotropic Compression Local Straight Edge Model 
Characteristics 


Hole Dia. (in.)/ 
Layup 

Nodes 

Degrees of 
Freedom 

Elements 

Local/Total 
Sole .on Time 
(CPU seconds) 

1.0/ 

[ ± 1 5jn]s 

3,239 

9,717 

576 

501.0/940.9 

3.0/ 

[±15 Jn ]s 

2,977 

8,931 

528 

488.1/925.8 

6.0/ 

[±15 Jn ] s 
[ ± 45 2 „]s 
[± 15 n /±45„] s 
[0„/±45„/90„]s 

2,715 

8,145 

480 

344.6/648.6 

342.0/693.0 

344.2/653.5 

347.0/776.2 


A concern associated with modeling curved geometries with finite elements is to have a sufficient 
number of elements in that region to reduce the amount of in-plane curvature modeled by the in- 
dividual elements. Through the selection of a greater number of spokes for the larger hole diam- 
eters, and thus more nodes, the in-plane curvature of the elements was kept low. Although more 
nodes were required in-plane around the local curved edge region for the larger hole sizes, there 
were fewer nodes required for the entire global model. The undeformed local curved edge models 
for the three different hole diameters are shown in Figure 64a), Figure 65a), and Figure 66a), while 
Table 9 contains quantitative information about the local models. Element arc lengths at a con- 
stant radial distance from the hole were varied to generate radial nodal lines from the hole edge to 
the global/local boundary around the hole. As with the cross-ply laminate models, the new lo- 
cations of the mid- side nodes for the elements were carefully selected to be exactly at the mid-point 
between the comer nodes. The resulting in-plane nodal distribution did not have the same angular 
distance between the individual nodes around the hole. 
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a) 




b) 


Figure 64. Curved Edge Models for |± IS^L with Six-Inch Diameter Hole: a) Undeformed b) 

Deformed 
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a) 




b) 


Figure 65, Curved Edge Models for |± ISjJg with Three-Inch Diameter Hole: a) Undeformed b) 

Deformed 
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Table 9. Three-Dimensional Angle-Ply and Quasi-Isotropic Compression Local Curved Edge Models 


Hole Dia. (in.)/ 
Layup 

Nodes 

Degrees of 
Freedom 

Elements 

Local/Total 
Solution Time 
(CPU seconds) 

1.0/ 

[±15 2n ] s 

6,432 

19,296 

640 

1,115.7/1,555.6 

3.0/ 

[±15 2n ] s 

8,040 

24,120 

800 

1,571.3/2,009.0 

6.0/ 

[±15 2n ] s 
[±45 2n ] s 
[± 15 n /±45 n ] s 
[0 n /±45 n /90 n ] s 

9,648 

28,944 

1,920 

2,330.3/2,634.3 

2,071.4/2,422.2 

2,090.3/2,399.6 

1,984.5/2,413.7 


5.3 Results and Discussion 


The two-dimensional global model was analyzed, and selected in-plane displacements were applied 
to the three-dimensional local models along the global/local interface. After determining the local 
nodal displacements, the strains and stresses at the integration points were calculated. Examination 
of solution times presented in Table 7 indicates that the two-dimensional to three-dimensional 
global/local finite element analyses presented were accomplished with a moderate amount of com- 
puter resources. 

The processor used to determine the displacement fields for all of these models was an iterative 
solver. The amount of computer resources needed by iterative solvers is often substantially less 
than those which contain conventional solution techniques. The time required for a solution with 
iterative techniques is dependent on the iterative method selected, the initial guess, and the condi- 
tion of the coefficient matrix, which is the system stiffness matrix. Since each two-dimensional 
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global model had a different stiffness matrix, models with the same total number of degrees of 
freedom required different amounts of computer time (Table 6 on page 101 ). No attempt was 
made to minimize the time to solve these systems, and the first set of parameters input to processor 
TTER" which yielded a solution with the correct accuracy level was used. The complete three- 
dimensional model for the [ ± 1 5 2n ] s laminate with a six-inch diameter hole did not converge to an 
admissible solution during repeated attempts to solve the system. Examination of the reported 
displacement field revealed that the elements were displaced in a manner inconsistent with the ap- 
plied load. It is believed that the stiffness matrix was strongly ill-conditioned, so that the tests for 
convergence incorrectly reported that the model was converged. With no available complete 
three-dimensional finite element solution, validation of the global/local modeling technique was 
completed, on a qualitative basis, by comparison to known information about the interlaminar 
stress fields in composite laminates. 

Angle-ply laminates whose stacking sequence is of the form [±©] s generate interlaminar stress 
fields through a different mechanism than cross-ply laminates. There is no Poisson's ratio mis- 
match between these laminae, while there is shear coupling at the interface. The strength and di- 
rection of the shear coupling is dependent on the magnitude and sign of the orientation angle of the 
plies. Unconstrained off-axis laminae with different ply angles tend to shear in different directions 
under uniaxial loading, but are constrained from doing so when in a laminate. Compatibility is 
enforced through the development of an interlaminar shear (r^) at the interface between the 
laminae. has been shown previously to be dominant and seemingly without bound at this lo- 
cation (Salamon, 1980). Previous studies of the interlaminar stress distributions for [±©] s lami- 
nates concentrate on this component (Pipes and Pagano, 1970), (Pagano and Pipes, 1971), (Dana, 
1974), (Pagano, 1974), (Hsu and Herakovich, 1977), (Rybicki and Schmueser, 1978), (Salamon, 
1980), (Bums et al, 1985), (Kassapoglou and Lagace, 1987). Most of these have focused on stresses 
at the straight free edge, but some did consider the interlaminar stress state around holes. 

More complex laminates formed with off-axis plies could take the form of a [±®i/±0 2 ]s stacking 
sequence, for example, or may be quasi-isotropic. Because these laminates include laminae of dif- 
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ferent angles, there may be Poisson's ratio mismatch or shear coupling or both at the ply interfaces, 
resulting in complex interlaminar stress fields to enforce compatibility and equilibrium through the 
laminate. However, certain general statements can be made about the character of these stress 
fields. The interlaminar normal stress was found previously to be maximum along the mid-plane 
for the [± 15 n /±45 n ] s and [0 n /±45 n /90 n ] s laminates. The interlaminar shear stresses (t^) is maxi- 
mum at the interface nearest the mid- plane for these stacking sequences, while the interlaminar 
shear stress ( t ^) exhibits maxima or minima at interfaces where it dominates. In addition, because 
the interlaminar stress field distribution is truly three-dimensional in character, the in-plane dis- 
tributions of each of these stresses change shape as well as magnitude for each ply interface. Such 
trends are reflected successfully in the interlaminar stress fields generated with the presented two- 
dimensional to three-dimensional global/local analysis method. 


5.3.1 Two-Dimensional Global Model Results 


The deformed grids for each hole diameter are shown in Figure 50b) through Figure 52b) for the 
[d=152n]s laminate. The deformed [±45 2n ]s plate with a six-inch diameter hole is found in 
Figure 53b), while the [± 15 n /±45 n ]$ and the [0 n / ± 45 n /90 n ]s deformed global models are in 
Figure 54b) and Figure 55b). All in-plane displacements were scaled by a factor of 10. Note that 
the vertical ends of each plate remained straight and their original length as specified by the global 
boundary conditions. The circular holes became noncircular, and there is a bulging of the previ- 
ously straight edges of the plates. The line of nodes two laminate thicknesses from the maximum 
global y-dimension is the global/local interface for the straight edge models, as it was for the 
cross-ply straight edge models. This boundary is curved outward from the hole region, and this 
curvature is transferred from the global model to the local through the applied displacement con- 
ditions. The square of nodes around each hole exhibits concavity on its vertical sides and convexity 
on its horizontal sides. These are the global/local interfaces for the local curved free edge model, 
and this deformation state will be applied to the local model. 


Angle-ply and Quasi-isotropic Laminates Under Compression 


125 



For a uniform end displacement, the reduced in -plane stiffness of those laminates with larger holes 
is reflected in the global deformation field, where the influence of the hole is experienced over a 
greater portion of the plate. Because the in-plane stiffnesses are also dependent on stacking se- 
quence, each lamination sequence will have a different in-plane displacement field for the same 
end- shortening displacement. Since the [± 15 2n ] s laminate has a higher x-direction in-plane stiffness 
than the [±45 2n ] s laminate, the same end-displacement deformed the [± 15 2n ] s plate less than the 
[ ± 45 2n ] s . Because the stiffness of the [ ± 1 5„/ ± 45„] s laminate is between these other laminates, the 
same end-displacement will cause more deformation than the [± 15 2n ]s laminate and less than the 
[±45 2n ] s laminate. By examining the global deformed grids in Figure 50b) through Figure 54, this 
is seen to be the case. No direct comparison between the interlaminar stress states are made because 
of the change in the in-plane laminate stiffness with stacking sequence. 

5.3.2 Three-Dimensional Local Model Results 


The straight free edge results shall be presented prior to those of the curved free edge. The success 
of the global/local analysis method shall be determined by comparing the trends of the interlaminar 
stress fields determined with the global/local finite element procedure to those from previous 
studies, since the complete three-dimensional model could not be solved. Preceding the straight 
edge model interlaminar stress results, two different methods of determining the straight edge dis- 
placement fields are presented and compared for the [db 1 5 2n )s plate with a six-inch diameter hole. 
One model was formed by selecting a smaller region of interest in the global model than the other, 
and it is this reduced model that will be compared to a full model for the straight edge region. Once 
it is established which model is preferable and why, the interlaminar stress fields at the straight edges 
near the different diameter holes will be considered. The effect of decreasing hole diameter in a 
finite-width plate for the [± 15 2 „] s lamination sequence is examined. The straight edge interlaminar 
stress fields for a plate with a six-inch hole and for the [ ± 4520]$, [ ± 1 5„/ ± 45„]s, and 
[0„/ ± 45„/90„]s laminates, respectively, are then studied. Since the region of interest for the straight 
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free edge is that nearest the hole, all of the presented information for the local models will be limited 
to that region. For this particular study that region is -3.5 < x < 3.5 in., where x is the local coor- 
dinate dimension. The in-plane displacements in the local straight edge models were scaled by a 
factor of 10, and the out-of-plane displacements by 100. In addition, their z-direction dimensions 
were multiplied by 5 to allow clearer viewing of the through-the-thickness effects. The curved free 
edge results will follow the same format: verification of the global/local model results, determi- 
nation of the influence of different hole diameters for the [± 15 2n ]s lamination sequence, and then 
the effect of different laminate orientations on the interlaminar stress field. While the entire region 
around the hole was modeled, only information at those angular positions between 0° <0 < 180° 
are presented since the data will be cyclic over the entire hole boundary. These are the integration 
points nearest the straight edge under examination and closest to the hole edge. 

5. 3. 2.1 Verification of Global/ Local Modeling Technique for the Straight Edge 

The primary feature of global/local finite element analysis methods is the option of selecting only 
those areas of interest for further discretization. With the two-dimensional to three-dimensional 
finite element analysis technique under consideration, the selection of this region and the displace- 
ment application procedure is not as straightforward as some analysts would prefer. Undeformed 
and deformed models for the two [± 15 2n ]s local straight edge models near a six-inch diameter hole, 
reduced and full, are in Figure 60a) and b) and Figure 61a) and b). Visual comparison of these 
deformed models suggests that there is a marked difference in the displacement of the free edges. 
Focusing in Figure 67b) on the displacements along the local x=3.5 in. plane in the full model, 
and on the ends of the reduced model, dramatically reveals the error in using a reduced model for 
this laminate and loading. The application of constant x-direction displacements along this model 
boundary inhibited shearing in the x-z plane, which is the dominant interlaminar effect for laminates 
with only off-axis plies. Because the global model cannot predict the exact through-the-thickness 
nature of the displacements and the application of constant displacements in the x-z plane is clearly 
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incorrect, local models of laminates with off-axis plies should not penetrate the free edge boundary 
layer. Furthermore, this will impact global/local stress analysis techniques for any application 
where the interlaminar stress state includes either x-z or y-z shearing: the local model may en- 
compass, but not partition, any region where shearing will be present. As a result, the local models 
along the straight free edge for any laminate with off-axis plies must be more extensive than those 
of cross-ply laminates under compression or tension. Further discussion of the straight edge stress 
analyses for laminates with off-axis plies will focus on the full local model, since the reduced local 
model displacement field is incorrect. 

Several attempts were made to verify the local straight edge stress fields with those found for a 
complete three-dimensional model for the plate with a six-inch diameter hole and a [± 15 2n ]s ori- 
entation. To allow comparison of the stresses at the same locations, the complete model had the 
same level of discretization as the local models. Attempts to solve the system with the Testbed 
were numerous, computationally intensive, lengthy, and unsuccessful. However, the degree of effort 
required to attempt this analysis proves the point that with a careful global/local analysis solution 
of otherwise intractable problems is a reality. Over 14 million coefficients form the complete 
three-dimensional system matrix with 4,608 elements, while the total number of nine-node plate 
elements for the global model was 576 and a total of 2,410 20-node brick elements were used for 
two local models. This represents a substantial saving of computational resources with this 
global/local technique over a complete three-dimensional finite element analysis. The validity of 
the global/local analysis is determined by comparing to previous studies in the open literature. 
Since finite element analyses using classical elements do not satisfy the traction-free conditions ex- 
actly, this precludes attempting to compare to stress field values generated from other methods. 
While there has been work which used finite elements in the past, close comparison to their results 
was not possible because of modeling differences between these and the current work. Since the 
global/local method and not the use of the finite element method is under scrutiny, until a complete 
three-dimensional analysis with the same discretization is completed, only the trends of the stress 
field can be validated. 
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Figure 67. Displacement Field Comparision Between Full and Reduced Models for |± 15^ near a 
Six-Inch Diameter Hole 
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Previous investigations have concentrated on the straight free edge of a [zt45 2n ] s laminate without 
holes, under uniaxial tension. Pipes and Pagano (1970) used a reduced exact elasticity formulation 
with a finite difference solution for a [±45 2n ] s laminate under uniaxial tension. The through-the- 
thickness distribution of the interlaminar shear stress (t xz ) was reported to approach a singular value 
at the interface between the layers. Values for the other interlaminar stress components {o z and 
T yz ) were presented along that interface with noticeably smaller magnitudes. The character of the 
stress distributions determined in their investigation is similar to that found with the current study. 
Wang and Crossman (1977) presented some results found with the finite element method, but not 
all of the trends between their study and the current one match exactly. However, some of their 
results do not match those of other researchers. Using a perturbation approach to solve the exact 
elasticity solution, Hsu and Herakovich (1977) found similar distributions to those in the current 
study. Rybicki and Schmueser (1978) modeled the ±6 layers as a single orthotropic layer, which 
eliminated the possibility of comparing stresses at the interface. Kassapoglou and Lagace (1987) 
presented a closed form solution for interlaminar stresses at a straight free edge which enforced a 
zero interlaminar normal stress at the interface. Both the in-plane and the through-the-thickness 
distributions at the interface were similar to those found by others, including this work. The 
global/local analysis technique determines the proper character of the interlaminar stresses along the 
straight edge for a particular x-location in the local model, and can be considered valid, 

53.2.2 Straight Free Edge Model Results 

The undeformed and deformed straight edge models for the [± 15 2n ] s plate with a six-inch hole are 
shown in Figure 60a) and b). From the deformed view of the model, the expansion of elements 
away from the laminate mid-plane along the straight free edge is seen to be confined within one 
laminate thickness of the free edge. Bowing of the elements nearest the hole away from the center 
of the plate is predicted, illustrating the influence of the nearby hole on the in-plane and out-of- 
plane displacement field. Careful examination of the deformed view confirms that the elements in 
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the center cf the local model a! > are expanding less through the laminate thickness. It is expected 
that these changes in the deformation field due to the hole will be reflected in the interlaminar stress 
distributions. The through-the-thickness distributions for the interlaminar normal stress (a z ) and 
interlaminar shear stress (r^) for this laminate are presented in Figure 08 a) and b). From 
Figure 68 a), the maximum interlaminar normal stress is near the laminate mid-plane, but com- 
paratively low in magnitude. The maximum stress for Ty* is found at the z = 0.025 in. interface in 
Figure 68 b) which is the interface between the 15° and -15° layers. Both of these figures underscore 
the fact that for a [± 15 2n ]s laminate these do not warrant further study because of their low mag- 
nitudes. 

The through-the-thickness and in-plane distributions at the interface for the dominant interlaminar 
stress component (t^) are presented in Figure 69a) and b) along the straight free edge in the 
[± 152 n]s plate with nearby a six-inch diameter hole. This stress component increases in value from 
zero at the mid-plane to a maximum at the ply interface, and decays to zero at the top free surface. 
Its character suggests that the interfacial values could continue to increase with further refinement 
of the model, indicating the possible presence of a stress singularity. The nearby six-inch diameter 
hole lowers the peak magnitude of this stress at the ply interface for a given x- location, but does 
not affect its in-plane character. Along the free edge, the lowest magnitude of is predicted to 
be at the centerline of the hole. As the distance from the straight free edge to the hole increases, 
and thus the effect of the hole is lessened, the the stress reduction becomes less prevalent. The 
expected limitation of the boundary layer to one laminate thickness in distance from the straight free 
edge is shown in the in-plane stress distribution in Figure 69b). 

General comments about the deformation and stress stress states at the straight edge for plates with 
three- and one-inch diameter holes follow, while the undeformed and deformed local models are 
presented in Figure 62a) and b), and Figure 63a) and b). In each of the deformed views of the local 
models, a thickening along those rows of elements nearest the straight free edge can be seen, with 
the greatest expansion across the final row of elements. Visual comparison of the undeformed and 
deformed meshes for all the hole diameters in the central area of the models reveals the diminished 


Angle-ply and Quasi-isotropic Laminates Under Compression 


131 





Figure 68. Through-the-Thickness Stresses for Local Straight Edge Model for |± IS^ls ncar Six-Inch 
Diameter Hole: a) c z b) 
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Figure 69. for Local Straight Edge Model for near Six-Inch Diameter Hole: a) 

Through -the -Thickness b) At the z* 0.025 in. Interface 
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effect of the hole on the deformation state for the smaller holes. The straight edge model near the 
one-inch diameter hole exhibits no in-plane curvature of the elements in Figure 63b), while those 
near the six-inch diameter hole in Figure 61b) are considerably bowed in the central portion of the 
model. 

Because of the similarities in the interlaminar stress distributions for the local models near different 
diameter holes, the individual results are not presented, but a direct comparison of the results fol- 
lows. The nearby hole affects all of the interlaminar stress distributions along the straight free edge, 
but its influence is clearest for the interlaminar shear stress (txz) at the z = 0.025 in. interface. The 
in-plane distributions of at the interface for the straight edge near the three- and one-inch di- 
ameter holes are similar to those near the six-inch diameter hole, except for the integration points 
nearest the free edge. The stresses over the localized region of interest for all three hole diameters 
are presented in Figure 70. Over this region, no hole effect is discernible for the plate with a one- 
inch diameter hole; as far as the straight edge interlaminar stress state is concerned, this plate does 
not have a hole. The three-inch diameter hole has some influence on the straight free edge 
interlaminar stress field, but not as marked as the six-inch diameter hole. Both of these holes de- 
press the stress magnitudes symmetrically in the central region of the plate. The straight edge values 
near the six-inch diameter hole are markedly depressed in this region; at the centerline, Tx* is nearly 
30% and 25% lower than those values for the straight edges near the one- and three-inch diameter 
holes, respectively. Near the ends of this region, the plate with a six-inch diameter hole does not 
recover the same stress as the other plates. This is attributable to the fact that the same end- 
displacement was applied to the plates, regardless of hole size. Since the in-plane stiffness for a plate 
with the six-inch diameter hole is lower, its in-plane stresses will be lower. Since the transition from 
the internal stress state to the traction-free edge state is over the same distance for these laminates, 
lower in-plane stresses will require lower stress gradients, and therefore, lower interlaminar stresses 
arise. The changes in the distribution and magnitude of the interlaminar stresses increase with in- 
creasing hole size, which suggests that the presence of the hole near the free edge reduces the 
intralaminar and interlaminar stress fields along the straight free edge. 
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Extensive examination of the influence of the in-plane compliances of laminates on their respective 
interlaminar stress fields have been conducted to determine which laminates develop large 
interlaminar stresses. The intention of this study was not to do further work in this area, but rather 
to demonstrate that the global/local modeling technique is able to determine these stress fields for 
a general lamination sequence. Two [±©]s laminates were considered prior to one with both of 
these angle of orientations: [±©i/±0 2 ] s . The first of these was a [± 15 2 „] s stacking sequence and 
the second, [±45 2n ] s , is considered next. 

When the global deformation fields were cursorily examined earlier, it was determined that the 
[±45 2n ]s laminate exhibited greater in-plane deformations than the [± 15 2n ]s laminate for the same 
applied end-shortening displacement. As seen when comparing the straight edge models in 
Figure 61b) and Figure 71b), the local model for the [±45 2n ] s laminate had a greater degree of 
deformation than the [± 1 5 2n ]s laminate, both in-plane and out. At the center of the plate the 
outward in-plane bowing of the elements is greater, and there is more element expansion through 
the laminate thickness for the [±45 2n ] s stacking sequence than the [± 15 2n ] s laminate. Near the 
center of the straight edge model, however, there is less thickening, as the six-inch diameter hole 
relieves the intralaminar stress and therefore, the interlaminar strain and stress state at the edge of 
the laminate. r„ is the dominant interlaminar stress component and is maximum at the interface 
between the plies, since this laminate is member of the [±0]s family. 

Unlike the [± 15 2n ] s laminate, the other interlaminar stress components are larger in magnitude for 
the [±45 2n ]s plate. Examination of the through-the-thickness distribution of a z in Figure 72a) 
reveals that it is nonzero and maximum along the mid-plane of the laminate. Its reported peak 
values are small, however, less than 5% of the normalization stress. The interlaminar shear stress 
(Tyj) is also small, reaching maximum absolute magnitude of 6% of the normalization stress, shown 
in Figure 72b). At the interface between the ±45° layers, exhibits opposing gradients. Aver- 
aging these integration point stress values to determine the interfacial values yields a near zero stress 
at the interface. Both of these stress components have variations in their magnitudes over this re- 
gion that are suggestive of a hole effect along the straight edge. 
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Figure 72. Through-the-Thickness Stresses at Straight Edge for |± 45^ laminate near Six-Inch Di- 
ametcr Hole: a) <r z b) 
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Figure 73a) is the through-the-thickness v * i laminar shear stress (t„) distribution for the 
[d=4 5 2n ]s laminate with a six-inch diameter he? j, while the in-plane t xz distribution at the interface 
between the ±45° layers for this laminate is presented in Figure 73b). One of the most noticeable 
features of the distributions in both of these figures is the symmetric decrease in magnitude of the 
stress about the centerline of the hole, which reaches a minimum magnitude at the centerline of the 
hole. This reduction is stronger than that found for the [± 15 2n ]s laminate and is attributable to 
the higher compliance of the [±45 2n ] s laminate. For each location in the localized region, the 
through-the-thickness distribution of for this laminate is similar to that for the previous stacking 
sequence examined, but with greater extremes in value. The in-plane distribution of at the 
z= 0.025 in. interface in Figure 73b) reflects the same information: its character is similar to the 
[± 1 5 2n ]s laminate with a larger peak magnitude for the [±45 2n ]s laminate, and the stress relief ef- 
fect due to the hole is stronger. These stresses vary from 0.5 times the normal in-plane stress of 
classical lamination theory to 0.15 at the center of the model. 

A laminate composed of ±15° and ±45° ply groups would be expected to exhibit characteristics 
of both of its ply groups. The undeformed and deformed models for the [± 15 n /±45 n ]s laminate 
are found in Figure 74a) and b). The through-the-thickness enlargement of the elements nearest 
the free edge appears less severe than the [± 15 2ri ] s laminate, but more severe than the [±45 2n ] s 
laminate. Bowing of the elements along the straight edge away from the center of the hole is clearly 
discernible for this lamination sequence, but is not as severe as that seen in Figure 71b) for the 
[±45 2n ] s laminate, and more than for the [±15 2n ] s laminate in Figure 64b). With higher dis- 
placement scaling factors than used here, deformations between the ply groups could be seen in the 
x-z and y-z planes. Because there is both a Poisson's ratio mismatch between the groups of layers 
and shear coupling for each group and between each group, the deformation field is complex and 
all three interlaminar stress components are expected to be nonzero for at least part of the boundary 
layer region. 

Examination of the through-the-thickness stress distributions near the straight edge for each com- 
ponent indicates which interfacial locations merit detailed study. All three of the through-the- 
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thickness stress distributions are presented in Figure 75a), b) and c). From Figure 75a), the 
locations of interest for the the interlaminar normal stress were identified as the laminate mid-plane, 
and the interface between the ±45° layers. Because the Poisson's ratio mismatch between the 
± 15° and ±45° orientations is not as strong as that between 0° and 90° layers, this stress compo- 
nent is lower in magnitude than for the cross-plies. The largest stress magnitudes were along the 
mid-plane and reached a maximum of less than 5% of the normalization stress. Reductions in the 
peak magnitudes of the stress attributable to the presence of the nearby six-inch diameter hole were 
determined to be greatest at the laminate mid-plane and diminished quickly through the laminate. 

The presence of steps in the through-the-thickness a z distribution is unsettling. Because of the 
dominance of x-z plane shearing and the fact that a 20-node element can only adequately predict 
a linear stress distribution, more compliant elements or a greater number of elements might have 
eliminated these steps. Recent work reported by Zahlan (1990) raised the possibility that more than 
two 20-node elements through a material layer are required some loading conditions. The use of 
two 20-node elements per ply is not an inherent assumption of this global/local analysis technique. 
Since the main thrust of the present work was to develop and verify this technique for general 
laminates, these steps are noted as a limitation of the element selection and not of the technique 
itself. 

From Figure 75b), the interface of interest for Ty, was identified as the interface between the -15° 
and 45° layers. Throughout this figure, the stress field exhibited several sign reversals, but at all 
locations its magnitude was small, less than 2% of the normalization stress. The presence of the 
hole has a moderate effect on the through-the-thickness distribution of this component which will 
be discussed later. 

The remaining stress component to be considered is i n . Because this laminate is composed entirely 
of off-axis plies, it is expected that this component would be dominant, and it is. A near zero stress 
can be found near the mid-plane and top of the laminate in Figure 75c). Also in this figure it is 
seen that high stress gradients appear on either side of each interface through the laminate thickness 
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at z = 0.0125 in., z = 0.025 in., and z= 0.0375 in. At the z= 0.0375 in. interface (between the ± 15° 
layers), the reported peak stress reached 30% of the normalization stress, while at the z — 0.0125 in. 
interface (between the ±45° layers), it reached 4% of the normalization stress. The ± 15° layers 
have a higher in-plane stiffness than the ±45° layers and will therefore resist the same end- 
shortening displacement to a greater degree, resulting in larger interlaminar stresses. Because 
stronger shear coupling also exists between these layers, the sharper and stronger peak at this 
interface is not unexpected. In addition, the stress reduction because of the nearby hole is strongest 
at this interface and decreases in influence as the z-location of the interfaces approaches the laminate 
mid-plane. It can also be seen that the shear stress values at the interface between the material layer 
groups are negative, while at the interfaces inside the layer groups the stresses are tensile. These sign 
reversals are consistent with the change in sign of the localized in-plane coefficient of mutual in- 
fluence through the laminate thickness. 

Previous interlaminar free edge stress distributions for the cross-ply and angle-ply laminates were 
affected by the presence of the nearby hole in a symmetric fashion about the centerline of the hole. 
However, two of the three in-plane stress distributions presented have asymmetric distributions. 
These are shown in Figure 76a) and b), which are the in-plane distributions of the interlaminar 
normal stress at the integration points nearest the mid-plane and at the ±45° interface, respectively. 
As the z-location increases from the mid-plane, o z decreases in value from 2.3% of the in-plane 
normal stress predicted by classical lamination theory (Figure 76a) to 1.0% at the far edges of this 
localized region (Figure 76b). From the through-the-thickness distribution in Figure 75a), steps 
in the interlaminar normal stress were noted for this laminate, so these reported values are some- 
what suspect. The in-plane distributions for one of the interlaminar shear stress components (r^) 
is presented in Figure 76c) and can be seen to have a symmetric distribution about the centerline 
of the hole at this interfacial location, although there is a strong asymmetry predicted for this stress 
component at other locations through the laminate. A peak magnitude of nearly 1% of the nor- 
malization stress was predicted by the analyses. 
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Figure 75. Through-the-Thickness Stresses for Straight Edge Model for l ± 15 n / ± 45,^ near Six-Inch 
Diameter Hole: a) o z b) v C)T « 
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Both of these stress components are small in magnitude throughout this laminate, and do not arise 
in-plane until the final row of elements at the free edge. The asymmetric distributions of these 
stresses were unexpected and could possibly be due to the changing load path around the hole and 
how this manifests itself in transverse shear through the laminate, or possibly the difficulty in sat- 
isfying the traction-free boundary conditions with classical displacement-based finite elements. 
Further mesh refmement, or even the selection of higher-order elements, would be advisable in this 
region, before firm statements can be made about the origin of this effect. 

The interlaminar shear stress (t^) in-plane distributions at the three ply interfaces are presented in 
Figure 77a), b), and c) in order of increasing z-location. For each of these distributions, there are 
large gradients in the elements nearest the free edge, suggesting again the existence of a stress 
singularity. Comparing these three distributions, their shapes and magnitudes change through the 
laminate. At the ±45° and ± 15° interfaces (Figure 77a) and Figure 77c), respectively), the char- 
acter of these curves for a given x-location is similar to that for the [±45 2n ]s and [± 15 2r J s laminates 
investigated earlier. The stronger interlaminar effect for the ± 15° plies is clearly discernible, as is 
the greater influence the hole. Since shear coupling between the plies is dependent on the order and 
magnitude of the ply angles at the interface, the effect in a given ±0 group will be dependent on 
the value of 0. At z = 0.025 in., the plies above and below this interface are -15° and +45°, re- 
spectively. The shear coupling between these plies is opposite to the other interfaces, which is re- 
flected in Figure 77b). This is the only interface whose distribution is not symmetric about the 
centerline of the hole for this stress component, so it would appear that the change in the load path 
around the hole and through the laminate thickness appears to affect the -15°/45° interface more 
than the other interfaces. Since the interlaminar normal stress at the mid-plane is directly dependent 
on those effects at the z= 0.0125 in. interface, and partially on those effects at the z= 0.025 in. 
interface, its asymmetry could be attributable to the -15°/45° interface. 

The complex stress state of the [± 15 n /±45 n ] s laminate was examined in detail to demonstrate two 
issues. As this lamination sequence is a mixture of the other angle-ply lamination sequences under 
consideration, the resulting stress field should resemble the other stress fields but still contain the 
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Figure 76. Stress Distributions for |± 15 n / ± 45 L near Six-Inch Diameter Holes a) a z at Mid-plane 
b) <r z at z “0.01 25 in. c) at Z“ 0.0125 in. 
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Figure 77. x„ for Straight Edge Model for [± 15 n /±45_L near Six-Inch Diameter Hole: a) At 

2*0.0125 in. b) At z* 0.025 in. c) At Z-0.03/S in. 
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effect of the Poisson's ratio mismatch not found in simpler lamination sequences. Trust in the 
global/local method is then gained before consideration of quasi-isotropic laminates. Quasi- 
isotropic laminates are used often in structural components, because their in-plane response is 
nearly isotropic in nature. Like the [± 15 n /±45„] s stacking sequence, quasi-isotropic laminates 
exhibit both shear coupling and Poisson's effects through their thickness due to the different layers 
that comprise the laminate. It will be clearly demonstrated with this global/local analysis that the 
interlaminar stress fields can be fairly complex, and not necessarily negligible. 

The undeformed and deformed local straight edge models for the [0 n / ± 45 n /90„] s laminate with a 
six-inch diameter hole are presented in Figure 78a) and b). The deformed view of this local model 
hints at the complex interlaminar strain state that exists. Because of the proximity of the six-inch 
diameter hole to the straight edge, the elements in the central area of the model bow outward away 
from the hole. Along the free edge, those elements which form the 90° plies immediately above the 
laminate mid-plane are severely compacted, with less compaction in the central portion of the 
model near the hole. Those elements nearest the top free surface represent the 0° layer and are 
greatly contracted toward the laminate mid-plane. With the displacement scaling selected, the final 
row of elements is contracted sharply downward, giving the illusion of a beveled edge. 

Examination of the through-the-thickness interlaminar stress distributions for this lamination se- 
quence reveals the interfaces of interest for the interlaminar normal stress (<r z ) and the interlaminar 
shear stresses, and t«, in Figure 79a), b), and c), respectively. Considerable through-the- 
thickness variation for each stress component was predicted. The a z magnitudes are strongly neg- 
ative at the mid-plane and increase in value with z-location, reaching a slightly positive value at the 
0°/45° interface before becoming zero at the top of the laminate. The six-inch diameter hole at the 
center of the plate reduces the magnitude of this component symmetrically about the centerline of 
the hole throughout the laminate half-thickness, with its strongest effect at the laminate mid-plane. 
Locations of interest for a z are identified from the through-the-thickness distribution as the mid- 
plane and all the ply interfaces. 


Angle-ply and Quasi-isotropic Laminates Under Compression 


148 




As shown in Figure 79b), the magnitude of changes from zero at the laminate mid-plane to a 
positive value at the -45°/90° interface, to a negative value at the 45°/0° interface, and back to zero 
at the top of the laminate. All of the ply interfaces will be investigated further. The reduction in 
the magnitude of the stress component attributable to the presence of the hole appears sharpest at 
the -45790° interface at z = 0.0125 in. 

It would be expected that the interlaminar shear stress (t«), which arises because of shear coupling, 
would be maximum through the laminate at the ± 15° interface, where shear coupling is also a 
maximum. In Figure 79c), this is seen to be the case. At the other ply interfaces, there are peaks 
in the stress which are smaller in magnitude. This stress component is zero at the mid-plane and 
top of the laminate. Two smaller peaks exist at the -45790° interface and the 0°/45° interface which 
are opposite in character to those at the ±45° interface. Again, the effect of the nearby hole is to 
reduce the magnitude of the stresses symmetrically about the centerline of the hole, and this re- 
duction is strongest at the interfacial location of the maximum stress. All interfaces are to be ex- 
amined in further detail for this component. 

In Figure 80a), the mid-plane interlaminar normal stress distribution for the quasi-iso tropic lami- 
nate resembles that for the cross-ply laminates. At the edges of the localized area, this component 
is reported to be 22% of the normalization stress, while stress reduction because of the nearby hole 
causes the stress to drop to 8% of the normalization stress. At the mid-plane, this stress reduction 
is symmetric about the centerline of the hole. In Figure 80b), the distribution at z = 0.0125 in., 
which is the -45790° interface, is presented. The shape of the predicted character of o z at this lo- 
cation is similar to that at the mid-plane but reduced in magnitude, with the peak stress at this lo- 
cation slightly more than half that found at the mid-plane. Again the presence of the hole in the 
laminate alters the stress distribution symmetrically about the centerline of the hole, and causes this 
stress to drop from 12% to 5% of the normalization stress. Figure 80c) is the in-plane o z distrib- 
ution at the interface between the ±45° layers [z — 0.025 in.). At this interface, there is a reversal 
in the sign of the peak stress value, and a further reduction in magnitude compared to the previously 
mentioned interfaces. The positive value here would suggest expansion of the layers, however, the 
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Figure 79. Through-the-Thickness Stresses for [0 n /± -M'Jg near Six-Inch Diameter Hole: a) 

b ) T yz c ) ff xz 
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magnitude of the stresses here is nearly 25% of the mid-plane values. This suggests that while there 
are tensile values, the deformation state may not strongly reflect this expansion because of the 
strong element contraction found at the the mid-plane. The stress distribution for the final interface 
of interest, the 0°/45° interface at z = 0.0375 in., is presented in Figure 80d). The in-plane distrib- 
ution here is similar to the z = 0.025 in. interface with slight tensile stresses at the free edge. As was 
shown in the through-the-thickness o z distribution, the magnitude of the stresses is reduced and the 
effect of the nearby hole is very small at this interface. 

Figure 81a), b) and c) contain the interfacial in-plane distributions of the interlaminar shear stress 
(Ty,) for the [0„/ ± 45 n /90„]s laminate near a six-inch diameter hole at z = 0.0125 in., z = 0.025 in., 
and z = 0.0375 in., respectively. The decay in magnitude of the stresses through the laminate half- 
thickness is evident from a cursory examination of these distributions. In Figure 81a), the distrib- 
ution of Tyt at the -45°/90° interface is presented. The reduction in stress due to the nearby hole 
appears symmetric over the centerline of the hole, with a 55% reduction from the maximum value 
presented. Of all the interfaces under consideration, this one has the highest magnitude for this 
component. The interlaminar effect is confined to one laminate thickness in depth from the free 
edge with high stress gradients over this region. However, those elements closest to the free edge 
contain stress values which appear quite different in character. The row of elements closer to the 
global/local boundary has stress values which would appear to grow without bound, but those 
nearest the free edge have a moderate slope. The possibility remains that an insufficient number 
of elements through a particular ply was used for this loading situation and mesh refinement studies 
should be completed prior to any firm conclusions. 

The next interface under consideration is the ± 45° interface at z-location z = 0.0125 in. The stresses 
at this interface are presented in Figure 81b). Comparing this figure to the previous one reveals the 
change in character of this stress component through the laminate thickness. The previous interface 
under consideration had steeper stress gradients through the boundary layer which were also op- 
posite in sign over the row of elements closest to the free edge. The maximum values are lower, 
reaching a high of 0.7% times the normalization stress at the far edges of the localized region as 


Angle-ply and Quasi-isotropic Laminates Under Compression 


152 




Figure 80. a z for |0„/ ± 45-/90 L near Six-Inch Diameter Hole: a) At Mid-plane b) At z« 0.0 125 

in. c) At z » 0.025 in. d) At z - 0.0375 in. 
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opposed to 0.15%. Again the stress reduction in the central region of the plate due to the nearby 
hole is evident. While not visible from the viewing angle selected, this reduction is not symmetric 
about the centerline of the hole. 

Figure 81c) contains the in-plane stress distribution of at the 45°/0° interface. A zero stress is 
reported until the final row of elements at the free edge, which exhibit steep negative gradients. 
From the viewing angle selected, it can be seen that there is a stress reduction effect attributable to 
the nearby six-inch diameter hole. However, it is not visible that this reduction in stress is not 
symmetric about the hole centerline. The shape of this distribution is different from the other 
interfaces in this laminate. 

The final interlaminar straight edge stress component to be presented is r„. Figure 82a) contains 
the in-plane distribution for this component at z = 0.0125 in., or at the -45°/90° ply interface. While 
there is some shear coupling between the 90° and -45° plies, this stress is zero over much of the local 
model at this interface. The rapid growth of this component over the row of elements closest to 
the free edge would suggest that either further discretization or the use of elements which satisfy the 
traction-free boundary condition exactly would alter its reported shape and magnitude. At the far 
edges of the localized region, a maximum magnitude of 2% of the normalization stress is realized. 
The view of the hole effect is obscured by the choice of viewing angle, but careful examination of 
Figure 82a) suggests that there is some reduction in this stress component. 

At z = 0.025 in., or the ±45° interface, r a is much larger in magnitude than the previous interface 
and exhibits a shape not unlike the [ ± 452„]s laminate. Shear coupling is strongest between the 
±45° layers than the other layers in this laminate, and therefore is the interface where the largest 
magnitudes of this stress component occur. Nonzero stress values are found only within one lam- 
inate thickness of the free edge, showing that the global/local analysis technique used to compute 
these results determined correctly that the expected interlaminar influence is confined to this region. 
From this figure, the strong influence of the nearby hole is seen to reduce this stress component 
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Figure 81. x for |0 n / ± 45 |l /90 n l K near Six-Inch Diameter Hole: a) At z~ 0.0125 in. b) At z* 0.025 

in. c) At z-O.oi/5 in. 
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symmetrically about the centerline of the hole from a peak value of 17% to 8% of the normaliza- 
tion stress. 

The remaining interface to be considered is that between the 07-45° plies. The in-plane values for 
Txz at this interface are presented in Figure 82c). As for the first interface considered for this stress 
component, the magnitudes reported here are small. Like the 907-45° interface, is nonzero only 
over the final row of elements along the straight free edge. It is interesting to note that there is an 
asymmetric distribution about the centerline of the hole which is visible in this figure. 

Examination of the interlaminar stress fields along the straight edge for the quasi-isotropic lami- 
nation sequence suggested quite definitively that while its in-plane response is isotropic, there is a 
strong interlaminar response. The interlaminar normal stress exists through the laminate thickness 
as does both of the interlaminar shear stresses. It was interesting to note that asymmetric distrib- 
utions about the centerline of the hole were found for several locations for this laminate and all of 
the stress components. Further examination would be required before firm conclusions of the or- 
igin of their asymmetric nature could be made. 

5. 3.2.3 Verification of Globalj Local Modeling Technique for the Curved Edge 

An unsuccessful attempt was made to verify the local curved edge model stress fields with a com- 
plete three-dimensional model of the [± 152n]s plate with a six-inch diameter hole. The validity of 
the results generated by the global/local finite element analysis procedure is therefore determined 
by a comparison of the trends in the interlaminar stress fields found with the global/local analysis 
technique to those described in the literature which were performed with other methods. Dana 
(1974) was one of the first investigators to apply three-dimensional finite element analysis to the 
determination of interlaminar stress distributions around holes in composite laminates. Tensile 
loads were applied to thick and thin [T45]» laminates with holes. While the thickness ratios are 
not identical between this and the present work, similar interlaminar stress variations about holes 
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Figure 82. Interlaminar Stress Distributions for |0 ft / ± 45 a /90 n ^ near Six-Inch Diameter Hole: a) 

at z“ 0.0125 in. b) at z* 0.025 in. c) T n at z* 0.0375 in. 
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in thin laminates were predicted. It was found that while the far-field laminate stacking sequence 
influenced the global stress state, the local interlaminar stress state was affected also by the effective 
lamination orientation around the hole. The effective stacking sequence relative to the hole surface 
changes around the hole because the global lamination sequence will be intersected at different an- 
gles by the curved edge. Even for laminates possessing only one type of interlaminar effect for the 
global lamination sequence, either Poisson's ratio mismatch or shear coupling, all three interlaminar 
stresses existed around the hole. Stress distributions through the laminate thickness were compared 
qualitatively for those angular positions presented. Agreement was good, but not exceptional. 
Because Dana's work was completed during the early 1970's, its efforts were limited by the available 
computational power, and the element discretization both in-plane and out are considerably coarser 
than the present study. Dana used three linear elements through plies, which was found by 
Vidussoni (1988) to be less accurate than two quadratic elements. Zahlan (1990) concluded that 
several quadratic elements are required through a ply thickness to determine the stress ffeld, casting 
further shadows on these results. Bums et al (1985) also used three linear elements through a ply 
thickness to investigate a broader range of laminate orientations with holes. Comparison of these 
results to the present work is difficult because only thick laminates were considered. The fact that 
the global/local modeling technique for the cross-ply composite laminate was as accurate as could 
be expected, the belief is that any errors in the global/local angle-ply or quasi-isotropic laminate 
modeling will come about because of discretization problems, not the global/local analysis proce- 
dure since similar trends were found by other researchers. 

S.3.2.4 Curved Free Edge Model Results 

When the local curved free edge models were generated for the angle-ply and quasi-isotropic lami- 
nates, the entire region around the curved edge was modeled. To give the reader an understanding 
of the computational effort involved with these local models, initially undeformed and deformed 
views are presented for each of the [± 15 2n ] s laminates in Figure 64 through Figure 66. From these 
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figures it is obvious that these models obscure the view of the hole edge, so several of the curved 
free edge models were cut in half and only the angular positions between 0° and 180° are presented. 
Each of these local models had its z-direction dimension multiplied by a factor of five to clearly 
reveal the through-the-thickness displacement and stress contour information. Any of the presented 
z-direction displacements were scaled by a factor of 100, while the in-plane displacements were 
scaled by a factor of 10. 

Figure 64a) and b) contains the undeformed and deformed view of the entire curved edge model 
for a [± 15 2n ]s laminate with a six-inch diameter hole. As noted previously, this model is extensive, 
with a large number of nodes and elements. The deformed global model had a rectangle of nodes 
bordering the curved edge region which exhibited both convexity and concavity, as shown in 
Figure 50b). These nodes are the global/local interfaces for the curved edge angle-ply local models, 
and it is these boundary conditions were applied to the local model and caused the hole to become 
noncircular, as seen in Figure 64b) for the local model. While the global model could only predict 
the noncircularity of the hole, the thickening and compaction around the hole was predicted by the 
local three-dimensional finite element model and presented for the deformed half-model view in 
Figure 83. The redistribution of the global compressive load because of the interruption of the 
global load path was manifested in the three-dimensional displacement or strain field in the local 
model. Furthermore, examination of this view reveals that the out-of-plane deformation field 
around holes in this laminate varies with angular position. This change in deformation appears to 
exhibit no symmetry about 8 = 90° along the curved free edge and is irregular enough to be attrib- 
utable to more than just the redistribution of the load path around the hole. Near the local x-axis, 
there is a strong compaction in those elements nearest the mid-plane of the laminates, while near 
8 = 90° there is both thickening and compaction of the elements through the laminate thickness. 
The three-dimensional character of this deformation state is not straightforward, as it is both a 
function of the change in load path and the effective stacking sequence around the hole. This 
precludes a strong reliance on intuition about the expected behavior around holes for angle-ply 
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laminates. Besides, if the deformation state is considered to be relatively complex, then the 
interlaminar stress state cannot be expected to straightforward either. 

In Figure 84, all three interlaminar stress components through the laminate half-thickness are pre- 
sented between 0 = 0° and 180° at those integration points nearest the curved edge of the hole for 
the [ ± 1 5 2n ]s laminate with a six-inch diameter hole. The interlaminar normal stress is nonzero 
around the hole, as seen in Figure 84a). Although there is no global Poisson's ratio mismatch for 
a [± 15 2n ]s laminate, there is a change in the effective Poisson's ratio between the plies around the 
curved free edge of the hole. Those edges not locally parallel or perpendicular to the global axes 
will experience local stacking stacking sequences other than [±0] s . At 0 = 0° and 90°, the localized 
lamination sequence will be a member of the [±0]s family, and a t will be zero at the free edge 
through the thickness. At 0 = 45° and 135°, the localized lamination sequence will be a [T©] s 
laminate which also has no Poisson's ratio mismatch. Examination of Figure 84a) reveals that at 
these angular positions, the thro ugh- the-thickness distribution contains lines of zero stress. For 
other locations around the hole, o z is maximum at the mid-plane and decays to zero at the top free 
surface with an anti- symmetric distribution about 0 = 90°. 

Figure 85 is a view of the deformed half-model with superimposed o t stress contours. The change 
in the sign of the stress contours near 0 = 90° over a small area in this figure was also shown in the 
through-the-thickness distribution in the previous figure. This method of presenting the stress 
contours on the deformed view of the model offers further insight into the stress state that exists 
around the hole. Those angular positions where there is enlargement of the hole edge are where 
there are concentrations of positive interlaminar normal stress contours, while those locations where 
the hole edge is compacted are where there are many negative contours. At 0 = 0° and 180°, there 
is compression of the elements both due to the in-plane and interlaminar effects. At the angular 
positions near 0 = 90°, the in-plane load causes expansion of the hole edge, while the localized 
Poisson's ratio mismatch causes both enlargement and compression of the elements nearest the hole 
edge. Further understanding of the complex interlaminar stress state around the hole requires ex- 
amination of the other stress components. 
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Figure 83. Partial Curved Edge Models for I ± I Sint around Six-Inch Diameter Hole: a) 

Undeformed b) Deformed 
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The through-the-thickness distribution of the interlaminar shear stress (t 9j ) is presented in 
Figure 84b), while the stress contours for this stress component are found in Figure 86. The 
dominant feature of the through-the-thickness distribution was the highly concentrated stress peak 
found at 0 = 90°. This high concentration of stress is reflected in the stress contour distribution as 
well. Because of the changing load path around the hole, high stresses are found in the global x- 
direction at this position, and consequently, large interlaminar stresses are not unexpected in this 
area. These large stress values reflect the strong tendency for the plies to move past one another 
in the plane parallel to the z-axis in this region. Examination of the deformed half-model reveals 
that at the same positions where the elements are moving most relative to another in the deformed 
grid, the stress contours for t S; are highly concentrated. 

The remaining interlaminar shear stress r n through-the-thickness distribution is presented in 
Figure 84c), and its stress contours are shown in Figure 87 on the deformed half-model. This is 
the stress component which reflects the tendency of the hole diameter change through the laminate. 
From Figure 84c), the magnitudes of this stress component at this radial location are small, with 
a complex shape which varies greatly with angular position and height in the laminate. This radial 
location is inside the laminate, but close to the hole edge where equilibrium considerations require 
that this stress component reach a value of zero although displacement-based finite elements cannot 
enforce this condition exactly. Because of the proximity of these integration points to the hole edge, 
the predicted stresses should be smaller in magnitude. The change in sign and magnitude with an- 
gular position mimics that of the interlaminar normal stress. For those locations identified earlier 
where the effective lamination sequence is symmetric and balanced, t„ is zero through the local 
model thickness; at intermediate values of 6, maxima and minima exist at the interface between the 
plies. Examination of the deformed model with superimposed stress contours reveals that the 
concentrated stress contours are centered on the interface at angular positions where the local model 
is changing from expansion to compaction in the z-direction or from compaction to expansion. 
As the laminate changes its z-direction movement, those elements on either side of the interface 
would cause the hole diameter to either reduce or expand, if permitted to move freely. Since dis- 
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Figure 86. t tz Contours for Partial Curved Edge Model for |± around Six-Inch Diameter Hole 
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placement compatibility must be enforced, this interlaminar stress component arises to constrain 
this motion. 

The influence of the different diameter holes on the interlaminar stress field around the hole is 
considered for the [± 15 2 „]s laminate for a six-, three-, and one-inch diameter hole in a 12 in. wide 
plate. From the through-the-thickness stress distributions for the interlaminar stresses around the 
six-inch diameter hole in Figure 84, locations of interest for each of the interlaminar stresses were 
identified. Incidentally, the distributions for the smaller diameter holes are not presented individ- 
ually, but as compared to the other holes since their through-the-thickness distributions were sim- 
ilar. The peak location and its magnitude changed as a function of the hole diameter for the stress 
fields, as determined by this global/local technique. 

Maximum values of the interlaminar normal stress were were predicted to exist at the integration 
points nearest the laminate mid-plane. Examination of o z for these three laminates reveals some 
interesting effects associated with the change in hole diameter (Figure 88) at the laminate mid- 
plane. The six-inch diameter hole stress distribution for this component is generally lower in 
magnitude than for the smaller holes, indicating that some redistribution of the load has occurred 
or that the stiffness of the plate is greatly reduced by the presence of this large hole. The stresses 
for the three-inch diameter hole are coincident with those for the six-inch diameter hole at the fol- 
lowing angular locations: O° < -0 < -25°, 75° <0 < 105°, and 155 <0<18O . However, the peak 
magnitude is largest for the three-inch diameter hole and its exact location is shifted further away 
from 0 = 90° than the six-inch diameter hole. Because of the reduced net section of the plate with 
the six-inch diameter hole, its compliance is less than the plate with a three-diameter hole. How- 
ever, the three-inch diameter hole is near enough to the straight free edge that some of the stress 
field is relieved. In the previous chapter, the three-inch diameter hole was identified as being closer 
to a 'critical' hole diameter for the cross-ply stacking sequence and the geometric dimensions of the 
plate. It is also so identified for the angle-ply laminates. The distribution for the one-inch diameter 
hole does not resemble those for the other sized holes as closely . The location of the peak magni- 
tude was shifted further from 0 = 90° than the three- and six-inch diameter holes. Between 
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0° < 6 < 45° and 135° <8 < 180°, the stresses for the one-inch diameter hole did not track those for 
the larger holes. Since the one-inch diameter hole which has no effect on the straight edge dis- 
placement or stress field, then the distribution of stress could not be expected to be very similar to 
the distribution of stress for holes which interact with the straight free edge. There are, however, 
fewer elements in the one-inch diameter curved edge model, and there is a possibility that the local 
model nodal discretization could have been further refined to produce a fully converged solution. 
Further examination of more fmely discretized models would determine which of these are domi- 
nating the stresses presented here. 

The next component under consideration is the interlaminar shear stress (t 9j ). From the through- 
the-thickness distributions in Figure 84b), its location of interest in the laminate was identified as 
the z = 0.025 in. interface. t Sz stresses for all three hole diameters are presented in Figure 89a) and 
can be seen to be quite similar in character. The three-inch diameter hole had the largest stress at 
the centerline of the hole, which was roughly 10% larger than the peak found for the six-inch di- 
ameter hole, while the one-inch diameter hole had the lowest value, nearly 70% smaller than that 
for the three-inch diameter hole. With the six-inch diameter hole displacement field interacting with 
the nearby straight edge, and the one-inch diameter hole having no interaction at all, the three-inch 
diameter hole could possibly be close to a critical hole size for this plate. Again the one-inch di- 
ameter hole does not have the same sharp character in its distribution as the other holes, raising the 
possibility of a less than fully converged local model. 

The peak magnitudes for r n are ten times less than those for the previously considered interlaminar 
shear stress and on the same order as the interlaminar normal stress, as can be seen by comparing 
Figure 89a) and b) and and Figure 88. Figure 89b) contains the stress values found for all three 
local models at z = 0.025 in. From the through-the-thickness distribution for the remaining 
interlaminar stress (t r ) in Figure 84c), the stress distribution is seen to contain interesting vari- 
ations. Near 6 = 90°, all exhibit both a negative and a positive peak and passes through zero at 
8 - 90°. The location of the peak stresses near 8 = 90° is closest to 8 - 90° for the six-inch diameter 
hole, and shifts further away with decreasing hole diameter. The exact location and magnitude of 
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Figure 89. Interlaminar Shear Stress Distributions at the Curved Free Edge for Different Hole Di- 
ameters: a) t #i at z- 0.025 in. b) T rt at z* 0.025 in. 
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the apex of this curve for the one-inch diameter hole is difficult to determine exactly because of its 
rounded nature. With increasing hole diameter, the maximum stress in this region increases in 
magnitude. Considering the stresses near 0 = 60° next, the effect of the change in hole diameter is 
less clear. The six-inch diameter hole has the lowest magnitude in this region, while the three-inch 
diameter hole has the largest. If the angular position of this peak stress for the three-inch diameter 
hole is taken as a reference, then the six-inch hole has its peak nearer 9 = 90° with the one-inch di- 
ameter hole closer to 0 = 0°. 

For the [± 152n]s laminate, the influence of increasing hole diameter in a finite-width plate is similar 
to the cross-ply laminates. The three- and six-inch diameter holes had the same distribution shape, 
while the six-inch diameter hole has a lower peak stress magnitude, suggesting that stress relief be- 
cause of the presence of the nearby straight edge did occur. Deviations from the shape of the 
stresses for the other holes were found for the one -inch diameter hole. Because of its small size, it 
has no influence on the interlaminar stress fields at the free edge, and therefore no stress relief. 
However, firm conclusions about its exact nature should not be made until a further study with 
different levels of in-plane nodal discretization are completed. Such a study would be more trac- 
table with the global/local method than with a full three-dimensional finite element analysis. 

To further understand the interlaminar stress state in laminated plates, a [ ± 45 2n ] s laminate with a 
six-inch diameter hole was modeled with this two-dimensional to three-dimensional global/local 
analysis technique. Most of the previous investigations of the interlaminar stress states at the 
straight free edge or around holes have included at least one [±45 Jn ]s laminate, so a wide body of 
information about the predicted stress states are available. While the entire local model was used 
during the analysis, only the partial deformed model is presented in Figure 90b). The in-plane 
compressive global loading causes the plate to expand in thickness, but the elements near the hole 
edge at times expand and contract depending on the localized displacement state. Near 0 = 0° and 
180°, the elements within one laminate thickness of the hole edge undergo strong compaction to- 
ward the laminate mid-plane. As the elements under consideration are further from this region, the 
degree of compaction lessens until the elements begin to expand. Unlike the [ ± 1 52„]s lami 
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this change in thickness is smooth and gradual. The ply angles are originally perpendicular to each 
other, and while their local effective lamination sequence will change as a function of the angular 
position, their relative orientation to one another remains 90° apart. The interlaminar effects due 
to the Poisson's ratio mismatch and the shear coupling between the layers can be described as 
functions of sin(0) and cos(0), so that relative ply orientation along the interface will cause the sines 
to be evaluated as cosines for the different plies. Since sine and cosine are smooth functions, the 
resulting deformation state would be expected to reflect this character. 

While the deformation state found for this laminate is smooth, the interlaminar stress state is not 
necessarily so. The through-the-thickness distributions of all three interlaminar stresses are pre- 
sented in Figure 91a), b) and c). Examination of Figure 91a), reveals that the interlaminar normal 
stress has a complex distribution which is maximum at the laminate mid-plane, which is consistent 
with that found for the other laminates considered. Near 0 = 90° where the in-plane stresses are 
high and changing over small distances, there are nonzero stresses near the interface between the 
plies. This is not consistent with previous results and may be attributed to the difficulty of 
matching the traction-free boundary condition at the edge with classical displacement-based ele- 
ments at this location. Examination of the o t stress contours on the deformed half-model in 
Figure 92b) concisely reveals that the negative stresses are concentrated in the elements which are 
nearest the mid-plane, while the positive stresses are found predominantly in the elements modeling 
the 45° ply. At the 'edges' of the partial model near 0 = 0° and 180°, no stress contours are found 
farther away than one laminate thickness from the curved edge of the model. Demonstrating again 
that the global/local modeling procedure has correctly determined the expected extent of the 
interlaminar stress state around the hole. 

The through-the-thickness r ei stresses are shown in Figure 91b), while the stress contours around 
the hole are in Figure 93. This stress component exhibits smooth changes in magnitude and shape 
as a function of angular position around the curved free edge. The through-the-thickness distrib- 
ution of has its peak values near 0 = 90° where the in-plane stress components are also large, 
with smaller peaks near 0 = 0° or 180°. At 0 = 45° and 0= 135°, this stress is zero through the 
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Figure 91. Through-the-Thickness Stresses for Partial Curved Edge Model for | ± 45^^ near Six-Inch 
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laminate. These locations are where the effective stacking sequence becomes a cross-ply laminate, 
which will exhibit no shear coupling between the plies. In addition, the deformed half-model has 
the largest concentration of stress contours in the region where the deformations are the largest. 
As seen in Figure 93, those elements near 9 = 90° at the ply interface must shear in the tangential 
direction as the elements in the 45° ply expand in thickness while those in the -45° ply contract, 
which is also the location of the peak stress magnitudes. 

In Figure 94, the dominant feature is the strong difference in value for t„ for those integration 
points above and below the interface. Because of the high value of i Bt over the laminate, the 
tendency for the layers to slide past one another at the hole edge is large, and is greatest in the region 
near 90°. To restrain this motion, the hole must change in diameter through its thickness, which 
results in significant differences in r n on either side of the interface. Near the plate where the in- 
plane loading and interlaminar effects are in opposition to each other, the difference in stress values 
are opposite for the interface near those elements modeling the region of the plate where these re- 
inforce one another. The jump in value across the interface is largest near 6 = 90°, because the stress 
field in this region has greater stress values for the other components, and therefore, larger differ- 
ences for T n result. Near 0 = 45° and 135°, this stress component is zero through the laminate 
thickness which are also the locations where the effective stacking sequence is a cross-ply laminate 
and the other interlaminar shear stress (t^) is zero. Examination of the partial model with the stress 
contours indicates the same information and trends. 

When the straight edge stress states were examined for a [± 15„/±45 n ]s plate with a six-inch di- 
ameter hole, a complex interlaminar stress state was determined. This is attributable to the 
Poisson's ratio mismatch between the groups of oriented plies and shear coupling at the center of 
each group and between the groups of plies. All three of the interlaminar stress components were 
examined in some detail. Studies of the interlaminar stress state around the six-inch hole for such 
a lamination sequence is difficult to interpret. As for the other laminates studied, the interlaminar 
stress state around the hole is dependent on the global stacking sequence, the effective stacking se- 
quence around the hole, the global applied load, and the localized load path around the hole. Even 


Angle-ply and Quasi-isotropic Laminates Under Compression 


176 



-.200- A 


190= B 


. 198= C 


.298= 0 


398= E 


-100= F 


Figure 93. t #z Contours for Partial Curved Edge Model for [±45^^ around Six-Inch Diameter Hole 



Angle-ply and Quasi-isotropic Laminates Under Compression 


177 


200 - ft 
150= B 
.108= C 
3588= 0 
9580= E 
108* F 
. 156= G 
.268= M 


Figure 94. Contours for Partial Curved Edge Model for around Six-Inch Diameter Hole 



Angle-ply and Quasi-isotropic Laminates Under Compression 


178 


those laminates whose global orientation exhibited but one type of interlaminar effect had a fully 
coupled interlaminar stress state around the hole and all three interlaminar stresses became signif- 
icant. Obviously, if the global stacking sequence will cause several interlaminar mechanisms be- 
tween its layers, the interlaminar stress state around the curved free edge of the hole will be 
extremely complex and difficult to interpret. All three of the interlaminar stresses will be expected 
to exhibit intricate distributions. 

To gain an understanding of the behavior of the interlaminar stresses around the hole, first the local 
deformation fields were examined. Figure 95a) and b) are the undeformed and deformed partial 
models for this stacking sequence and hole diameter. When a global compressive loading is applied, 
the plate will expand through its thickness for a positive laminate through-the- thickness Poisson's 
ratio, which is the case for this laminate and all of the other laminates under consideration and is 
seen in this figure. As determined for several of the other stacking sequences, those elements nearest 
the hole edge at 0 = 0° and 180° will deform toward the laminate mid-plane. As 0 increases, this 
compaction lessens and reaches zero at 0 = 45°. When 45° < 0 < 135°, the elements nearest the hole 
edge can experience compaction and enlargement through the laminate thickness. Following the 
top row of elements nearest the hole edge in Figure 95b) reveals this effect for this layer of elements, 
while the deformation of the normals to the plies suggests that this deformation state is not constant 
through the laminate. 

The trends in the deformation field are reflected in the through-the-thickness distribution of the 
interlaminar normal stress, which is presented in Figure 96a). The through-the thickness distrib- 
utions for the interlaminar shear stresses (t* 2 and j n ) are found in Figure 96b) and c). As expected, 
all of the stress components are complex. The interlaminar normal stress has a complex variation 
through the thickness along the hole edge, because each angular position now experiences Poisson's 
ratio mismatches at three interfaces for those positions away from 0= 0°, 90° and 180°. Concen- 
tration of the stress at these different z-locations in the laminate is shown for the stress contours in 
Figure 97. At 0 = 0° and 180°, the through-the-thickness distributic :een to be similar to the 
straight edge for this laminate, while at 0 = 45° and 135°, there is no interlaminar normal stress 
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through the laminate. Examination of the stress contours superimposed on the deformed grid 
suggests that those areas with the greatest change in thickness are those with the largest stress gra- 
dients, which is not unexpected. At 0 = 0° and 180°, the largest concentration of negative stresses 
are found along the mid-plane in the region with the strongest element compaction. Near 0 = 90°, 
some element expansion was found as well as some positive stress values at the 15 /-45 interface. 
Examination of both the through-the-thickness stress distributions and stress contours for this stress 
component clearly indicates that the stress distribution is not symmetric about the hole centerline. 
A cursory examination of the stress contours of the [± 15 2n ]s 2 nd [±452n]s laminates show that the 
distribution for the [± 15 n / ± 45 n ] s laminate is made considerably more complex by the enforcement 
of compatibility between the 15°/-45° layers. 

Examination of the through-the-thickness interlaminar shear stress ( t Bz ) in Figure 96b) reveals that 
its distribution resembles more closely a combination of those for the [± 15 2n ]s and [±45 3n ] s lam- 
inates. The [± 1 5 2 n]s laminate had a highly concentrated stress near 0 = 90°, and this concentration 
can be found at the interface for the ± 15° layers of the [± 15 n /±45 n ]s plate, but shifted slightly 
toward 0 = 0°. The [±45 2n ]s laminate, on the other hand, had a more gradual stress distribution 
as a function of angular position, and this relative smoothness can be seen in the stress distribution 
for the [± 15 n /±45 n ]s laminate. At the 15°/-45° interface, this interlaminar shear stress is negative 
in value and shifted away from 9 = 90° toward 0 = 180°. While this is not clear from the through- 
the-thickness distribution, it can be seen in the stress contours on the deformed partial model in 
Figure 98. For the [±45 2n ]s laminate, there were additional stresses near 0 = 0° and 180°. This 
character was determined at the ±45° interface for the [± 15 n / ±45 n ]s laminate. 

For both of the previous angle-ply lamination sequences, the r„ distribution around the hole and 
through-the-thickness was complex in shape and difficult to describe. As expected, this stress 
component is no less so for the [± 15 n /dt45 n ]s laminate. Its through-the-thickness and stress 
contours around the hole are found in Figure 96c) and Figure 99, respectively. The through-the- 
thickness distributions resemble none of the distributions shown for the previous lamination se- 
quences considered and are difficult to interpret. As the effective stacking sequence changes around 
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Figure 96. Stresses for Partial Curved Edge Model for | ± 15 B / ± 45,^ around Six-Inch Diameter Hole 
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the hole, strains are developed at each interface to maintain compatibility of the layers. Each an- 
gular location which has concentrations of stress contours corresponds to a position which has 
considerable deformation along the row of elements along the top of the plate. However, the lo- 
cation of the stress concentrations through the laminate thickness changes with angular position. 
Several of the peaks found in the through-the-thickness distributions are seen to be at slightly dif- 
ferent angular positions and at different heights in the laminate. Most of the strong peaks were 
at the ± 0 interfaces, while the 15°/-45° interface experienced a lower stress magnitude over a larger 
region around the hole. While it is not clear from the deformed partial model stress contours, it 
appears that the distributions were somewhat anti -symmetric about 0=90°. The positive values 
were at the ±45° interface on one side of 0 = 90° and at the ± 15° interface on the other side of 
0 = 90°. Likewise, the negative shear values were found on either side of 0 = 90° for the different 
interfaces. 

The [±0,/±0 2 ] s lamination sequence study of interlaminar effects around holes is difficult to in- 
terpret because all of the three interlaminar stresses exist and have strong gradients around the hole 
at several interfaces. Each interface around the hole is affected by Poisson's ratio mismatch between 
the ply groups, and shear coupling within and between the ply groups for most of the curved edge. 
The resulting interlaminar stress distributions are difficult to interpret at best, and will be for any 
laminate with several ply interfaces through its thickness. 

Quasi-isotropic laminates are similar to the [ ± 1 5 n / ± 45 n ] s laminate since there are several ply 
interfaces where Poisson's ratio mismatch and shear coupling can possibly occur. For a 
[0 n /±45 n /90 n ]s plate with a six-inch diameter hole, each ply interface then becomes of interest. 
The undeformed and deformed partial models for this laminate are presented in Figure 100a) and 
b). Around the hole, thickening of the model due to the in-plane loading is gradual and smooth. 
However, the deformations in the boundary layer around the hole enforces compaction of the ele- 
ments for all angular positions which varies in intensity. Close examination of the deformed model 
reveals that at the laminate mid-plane all of the elements are compacted. In the ±45° layers, there 
is both expansion and contraction, and in the 0° layers, enlargement of the elements can be detected 
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around the hole edge. Sliding e elements relative to one another in the x-y plane occurs at 
several locations around and thr^^h the laminate. Definitive changes in the hole diameter through 
the laminate can be seen at the 0° and 90° locations because of the cut in the local model. 

In Figure 101, the through-the-thickness distributions for all three interlaminar stress components 
are presented. Figure 101a) is the interlaminar normal stress distribution, while the two 
interlaminar shear stresses, r yz and T nt are in Figure 101b) and c), respectively. From 
Figure 101a), the interlaminar normal stress can be seen to be large and mostly negative, except in 
the ±45° plies. Maximum magnitudes were found near the laminate mid-plane. Figure 102 is a 
view of the partial deformed model with a z contours superimposed. The positive stresses are con- 
fined to the ±45° plies and near 0 = 90°. However, the magnitudes of these positive stresses are 
small in comparison to the negative stresses found throughout the curved edge region. As a con- 
sequence, there is no expansion of the elements around the hole attributable to the interlaminar 
strain state. There is, however, a lessening of the compaction of elements around the hole. At 
0 = 0° and 180°, the elements are strongly compressed toward the laminate mid-plane, while the 
through-the-thickness stress distribution indicates significant stress values only in the 90° layers at 
these locations. 

Figure 101b) is the through-the-thickness t* z stress distribution. At the 0°/45° interface, the peak 
magnitudes are slightly larger than those for the ± 45° interface and considerably larger than at the 
-45°/90° interface. This change in value is attributable to the amount of load carried by the re- 
spective layers. For the same end-shortening, the 0° layers will experience a higher in-plane stress 
than the |45°| or the 90° layers, and so its interlaminar stress values are expected to have larger 
magnitudes. Furthermore, the location for all of these peaks is near 0 = 90°, which is a location of 
concentrated in-plane stress for even the global model. The deformed partial model with the stress 
contours in Figure 103 reveals the location of some of these peaks more clearly through the lami- 
nate. At 0 = 0° and 180°, the sign of the stresses is opposite in this figure. This shows that the 
sliding of elements is in the opposite direction. Near 0= 135°, the compaction of the elements is 
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Figure 100. Partial Curved Edge Model for [0 n / ± 45^/90,,^ around Six-Inch Diameter Hole: a) 

Undeformed b) Deformed 


Angle-ply and Quasi-isotropic Laminates Under Compression 


188 




Figure 100. 


uSaSr?) r " •"»* S-'-cl. Note 


a) 


Angle-ply and Quasi-isotropic Laminates Under Compression 


188 







. 680 = A 

B 

•300= C 
208= D 
180= E 
108= F 
260 = 6 


Figure 102. <r ; Contours for Partial Curved Edge Model for Quasi-isotropic Laminate around Six-Inch 

Diameter Hole 



Angle-ply and Quasi-isotropic Laminates Under Compression 


190 


strong and results in the shearing of elements in one direction. However, near 9 - 90°, the elements 
are experiencing less compaction and the sign of the stresses is opposite. 

Figure 101c) is the through-the-thickness stress distribution for the remaining interlaminar stress 
component, r re . This stress distribution is difficult to interpret. Its magnitudes are larger for this 
laminate than any of the other lamination sequences considered in this study. On either side of 
6 = 90°, the stresses peak and change direction in the ±45° layers. In the 90° layers, there are large 
values in those same locations where the element compaction is greatest. Obviously, those elements 
are shifting through the laminate and changing the hole diameter slightly. As can be seen in 
Figure 104, this areas are coincident with those where the thickening and contraction of the ele- 
ments through the laminate thickness are changing. 


5.4 Conclusions for In-Plane Compression of Angle-Ply 
and Quasi-Isotropic Laminates 


Several conclusions can be drawn from the analyses discussed, both about the state of stress around 
large holes in thin laminates and to the application of two-dimensional to three-dimensional 
global/local stress analysis in laminated composites. Considerable savings in computer time and 
storage can be achieved by using a global/local approach for analysis of laminated composites. It 
should be noted that use of the global/local method made the analysis of a thin laminated plate 
possible, as the full three-dimensional model would not converge. 

During this study several conclusions were drawn about interlaminar stress and strain states for 
angle-ply and quasi-isotropic laminates. Interlaminar stress fields for straight and curved free edges 
were studied for varying stacking sequences and hole size. 
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Interlaminar stress fields in laminates with off-axis plies are more complex and dependent on 
stacking sequence than cross-ply laminates, because of the wider range of possible ply angle com- 
binations, and therefore, interlaminar effects. Angle-ply laminates of the [± 0]s orientation do not 
experience a Poisson's ratio mismatch between the layers, but do have coupling between extension 
and shearing. The degree of shear coupling present, in turn, is dependent on the angle 0 in the 
laminate stacking sequence. 

For a common end-shortening, because the compliances for the laminates considered can vary 
greatly, different in-plane stress levels result, which influence the interlaminar stress state. The 
[± 152n]s laminates studied are less compliant than the [±45 2n ]s laminate. Interlaminar stress fields 
along straight edges in proximity to a hole arise because of shear coupling for [±©]s laminates. 
The [ ± 1 5 2n ] s lamination sequence has a stronger shear coupling effect than the [ ± 45 2n ] s laminates, 
which was reflected in the interlaminar stress fields. 

Combining ±0 layer groups together to form a [±0!/±0 2 ] s laminate immediately changes the 
interlaminar stress state to one resembling quasi-isotropic laminates. Between each group of layers, 
there is a Poisson's ratio mismatch, while shear coupling can occur at every layer interface. At a 
straight edge, the interlaminar normal stresses will be maximum at the mid-plane. At a nearby 
interface, the interlaminar shear stress will be maximum. The interfacial location of the maxi- 
mum r a is dependent on the exact stacking sequence considered. 

The effect of a hole on the nearby straight edge, if there was an effect, was to reduce magnitude of 
the interlaminar stresses for all laminates and hole diameters considered. The one-inch diameter 
hole appears to have no effect on the straight edge for all of the stacking sequences under consid- 
eration, while the six-inch diameter hole greatly influences the distribution character. The three- 
inch diameter hole stresses were determined to be a combination of those for the one- and six-inch 
diameter holes. 
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The interlaminar stress distribution around holes is more complex, and less intuitive. Like the 
stresses near the straight edge, it is dependent on the stacking sequence. However, as the angular 
position around the hole changes, so does the effective stiffness and stacking sequence. This affects 
the interlaminar stress fields greatly, and is why, although no interlaminar normal stress develops 
for [±0] s stacking sequences at the straight edge, it does around the hole. The local Poisson s 
ratio mismatch and shear coupling change in value, which affects the stress fields. Some of the 
lamination sequence distributions are less complex than others, but intuition can lead the analyst 
astray. 

The effect of the different hole sizes on the interlaminar distribution around holes is more involved 
and highly lamination sequence dependent. General comments are that the three- and six-inch di- 
ameter stress distributions resembled each other closely. The locations of stress peaks were gener- 
ally shifted slightly and not at the same location as for the cross-ply laminates. In portions of the 
plate farther from the straight edge, the distributions could be coincident, but nearer to it there were 
more deviations. The three-inch diameter hole exhibited larger peak values for most of the stress 
components than the six-inch diameter hole, and this would suggest that a critical hole size exists 
in the laminate somewhere near three inches. It is believed that conclusive statements about the 
effect of the one-inch diameter hole in the hole region are premature until the issue of whether the 
deviations in the stress distributions are because of the lack of stress relief, or merely a discretization 
deficiency. 
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6.0 Conclusions - Global/Local Analysis 


The two-dimensional to three-dimensional global/local finite element analysis technique presented 
was applied to several laminates with holes under compressive loading. Cross-ply, angle-ply and 
quasi-isotropic stacking sequences were modeled with three different hole diameters. Plate length- 
to-thickness ratio was 180, plate width- to -thickness was 120, and hole diameter- to-thickness values 
were 60, 30 and 10. Interlaminar stress fields at the straight edge near the hole and around the hole 
were examined in detail. Few previous investigations have involved these aspect ratio ranges, in 
part due to the prohibitive computational requirements of three-dimensional finite element analysis. 
Considerable saving in computer time and storage was achieved by using a global/local approach 
for analyses of laminated plates. The numerical results that were presented could have been more 
detailed and accurate through the use of more refined meshes, but the purpose of this study was to 
develop and verify the applicability of the procedure to the analysis of laminated composites 
structures. One of the primary goals of this global/local analysis method is to reduce the resource 
requirements of detailed analyses by performing several moderate analyses rather than one ex- 
tremely large analysis. 

For the cross-ply laminates, comparison between a complete three-dimensional model and local 
models of a plate for a specific hole diameter and stacking sequence was excellent along the straight 
free edge, and reasonable along the curved edge for all interlaminar stress components. Better 
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agreement is believed possible with an improved in-plane distribution of elements in certain lo- 
cations around the hole. The global/local analyses were able to determine the effects of clustered, 
dispersed, and reversed plies from the same in-plane displacement field on the three-dimensional 
interlaminar stress field. Conclusions about the hole effect in a finite plate on the interlaminar stress 
fields at a nearby straight edge and around the holes were made, and presented at the conclusion 
of that chapter. 

The global/local displacement application procedure for the straight edge cross-ply models proved 
to be a special case where prior knowledge of the strain field (y« = 0) allowed the free edge boundary 
layer to cross the boundary of the local model. All symmetry planes were exploited to reduce both 
the straight and curved edge models to their smallest possible form. Care must be taken in the 
displacement procedure to not inhibit boundary layer development with kinematically inadmissible 
displacements along the symmetry planes, and the 'cut' through the boundary layer. 

For laminates with off-axis plies, the three-dimensional nature of the strain and stress fields preclude 
the application of constant displacements in any direction through the boundary layer as was per- 
missible in one plane for the cross-ply laminates. Thus, the straight edge models had to extend from 
one end of the plate to the other. The error in attempting to reduce the size of the local angle-ply 
models by cutting a boundary layer was explored. In order to model the global load application 
conditions with the reduced local models, assumptions about the displacement field were made at 
the cuts. However, the only possible assumption was kinematically inadmissible and inaccurate 
displacement and strain fields developed. The local curved free edge models for the cross-ply and 
angle-ply laminates are similar in character. However, the cross-ply laminates under in-plane 
loading contain several symmetry planes which were exploited that do not exist for the other lam- 
inates or for more general loading conditions. Thus, the size of the angle-ply models around the 
hole are nearly four times that of the cross-ply laminates. 

An attempt was made to correlate the global/local analyses of an angle-ply laminate to a complete 
three-dimensional model. However, examination of the displacement fields from the analysis indi- 
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cate that this model had not converged. Trends from previous analyses indicate that the 
interlaminar stress fields at the straight free edge and around the hole are reasonable and at least as 
accurate as those analyses. The fact that the three-dimensional model was not as tractable as be- 
lieved is a sign that the global/local analysis technique can be applied to more intricate problems 
to assist in the analysis and design of structures. 

Cautions and general comments about the global/local analysis procedure should be made. Ap- 
plication of the global/local technique proved to be successful when some basic rules are followed. 
The global analysis must model the global response accurately. Generally, the global/local inter- 
faces should be located in regions with relatively uniform stress fields, making an extremely fine 
global grid unnecessary, and its exact placement was found not extremely critical as lonj it is 
outside the region of local influence. Of course, the global/local interfaces should be . to 

the local effect as possible to reduce the model size without overconstraining the devel of 

the three-dimensional displacement or stress fields. With proper mesh design, the mode] 5^. and 
computer resources can be kept low. A requirement on the local model is that it determines the 
local effect accurately. Certain short-cuts are permissible, but are dependent on the three- 
dimensional displacement fields. Through the global/local method, it is possible to tailor problem 
size to the accuracy required, and to more specifically model local effects in thin laminated plates. 
When applying such a technique, care must be taken by the analyst, to have an understanding of 
both the two-dimensional and three-dimensional finite element analyses limitations and capabilities, 
and the two-dimensional stress state, if not the three-dimensional stress state. Without such an 
understanding, an inappropriate assumption might be made and the global/local analysis would 
yield inappropriate displacement and stress fields. 

Global/local analyses can drastically improve the efficiency of computerized structural analysis and 
can provide sizable savings by circumventing the need to perform expensive analyses of non-critical 
regions. In addition, the design of the Testbed and the work of the Computational Mechanics 
Branch assisted the present effort tremendously by having specialists create modular proa sors 
which take advantage of local memory management schemes, utilize sparse matrix storage, and 
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contain efficient solution techniques. The result is a tool which minimizes the effort of the struc- 
tural analyst in determining structural response, thus permitting more complicated analyses, that 
were previously intractable, to be completed. 
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